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Abstract 

We settle the issue of well-posedness for the Dirichlet problem for a higher order elliptic 
system C(x,D x ) with complex- valued, bounded, measurable coefficients in a Lipschitz domain 
Q, with boundary data in Besov spaces. 

The main hypothesis under which our principal result is established is in the nature of best 
possible and requires that, at small scales, the mean oscillations of the unit normal to dfl and 
of the coefficients of the differential operator C(x, D x ) are not too large. 

1 Introduction 

A fundamental theme in the theory of partial differential equations, which has profound and in- 
triguing connections with many other subareas of analysis, is the well-posedness of various classes 
of boundary value problems under sharp smoothness assumptions on the boundary of the domain 
and on the coefficients of the corresponding differential operator. In this paper we initiate a pro- 
gram broadly aimed at extending the scope of the agenda set forth by Agmon, Douglis, Nirenberg 
and Solonnikov (cf. [4 J , [45J , [46J ) in connection with general elliptic boundary value problems on 
Sobolev-Besov scales, as to allow minimal smoothness assumptions (on the underlying domain and 
on the coefficients of the differential operator). Our main result is the solvability of the Dirichlet 
problem for general higher order elliptic systems in divergence form, with complex- valued, bounded, 
measurable coefficients in Lipschitz domains, and for boundary data in Besov spaces. In order to 
be more specific we need to introduce some notation. 

Let m, I £ N be two fixed integers and, for a bounded Lipschitz domain f2 in R™ (a formal 
definition is given in §6.1) with outward unit normal v = (y\, v n ) consider the Dirichlet problem 
for the operator 

C(X,D X )U:= D a {A aP {X)D^U) (1.1) 

\a\ = \f3\=m 

i.e., 



*2000 Math Subject Classification. Primary: 35G15, 35J55, 35J40 Secondary 35J67, 35E05, 46E39. 
Key words: higher order elliptic systems, Besov spaces, weighted Sobolev spaces, mean oscillations, BMO, VMO, 
Lipschitz domains, Dirichlet problem 

The work of authors was supported in part by from NSF DMS and FRG grants as well as from the Swedish National 
Science Research Council 



1 



D a (A a/3 (X)D^U) = for X ESI, 

\a\ = \/3\=m 

(1.2) 

— -r = 5fc on 90, < k < m - 1. 
- en/* 

Here and elsewhere, D a = (— id/dx\) ai ■ ■ ■ (—id/dx n ) a " if a = (ai, ...,a n ). The coefficients ^4 a/ 3 
are I x Z matrix-valued functions with complex entries satisfying 

Y IIA^IU^o) < «i (1-3) 

|a|=|/3|=m 

for some finite constant K\, and such that the coercivity condition 

K X E (^(^W(X),D Q W(X))dX>K 1] ||£> Q W||! 2(Q) (1.4) 

|a|=|/3|=m |a|=m 

with kq = const > holds for all C'-valued functions IA E C^°(S7). Throughout the paper, "Stz 
denotes the real part of z E C and (•, •) stands for the canonical inner product in C\ 

Since, generally speaking, v is merely bounded and measurable, care should be exercised when 
defining iterated normal derivatives. For the setting we have in mind it is natural to take d k /dv k := 
(Sj=i £jd/dxj) k \^ =v or, more precisely, 

^ := i* ^-u a Tr[D a U], < k < m - 1, (1.5) 

|a|=fc 

where Tr is the boundary trace operator and v a := v® 1 ■ ■ ■ v% n if a = (a±, a n ). With p(X) := 
dist (X, <9$7) and p £ (l,oo), a E (— 1/p, 1 — 1/p) fixed, a solution for 1)1.2)1 is sought in PF™' a (Sl), 
defined as the space of vector-valued functions for which 

[ \D a U{X)\ p p(X) pa dXy /P < oo. (1.6) 
0<|a|<m n 

In particular, as explained later on, the traces in 1)1.5(1 exist in the Besov space Bp(dQ), where 
s := 1 — a — 1/p E (0, 1), for any W E W^"' a (f2). Recall that, with da denoting the area element on 

dn, 



geB s p (dn)^\\g\\ B s {m) :=\\f\\ Lp{m) + ([ f 'ff Y{ £}1 da x da Y ) 1/P < oo. (1.7) 

The above definition takes advantage of the Lipschitz manifold structure of <9S7. On such manifolds, 
smoothness spaces of index s E (0, 1) can be defined in an intrinsic, invariant fashion by lifting 
their Euclidean counterparts onto the manifold itself via local charts. We shall, nonetheless, find 
it useful to consider higher order smoothness spaces on 9S7 in which case the above approach is no 
longer effective. An alternative point of view has been developed by H. Whitney in [S'2j where he 
considered what amounts to higher order Lipschitz spaces on arbitrary closed sets. A far-reaching 
extension of this circle of ideas pertaining to the full scale of Besov and Sobolev spaces on irregular 
subsets of W 1 can be found in the book [26J by A. Jonsson and H. Wallin. 
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For the purpose of this introduction we note that one possible description of these higher order 
Besov spaces on the boundary of a Lipschitz domain Q C M. n and for m£N,p£ (1, oo), s G (0, 1), 
reads 

B™- 1+s (dQ) = the closure of \{D a V|an)|a|<m-i : V G C °°(lR n )} in B s p {89) (1.8) 

(we shall often make no notational distinction between a Banach space X and 3L N = 3t © • • • © X 
for a finite, positive integer N). A formal definition along with other equivalent characterizations 
of B™~ l+s {dVL) can be found in §6.4. 

Given (|1.5jl - l|1.6|) . a necessary condition for the boundary data {<?fc}o<fc<m-i in (|l-2j) is that 

there exists / = {f a }\ a \<m-i € B™~ 1+s (dQ) such that 

/c! (1.9) 
9k = i k — r v a f a , for each < k < m — 1. 
a! 

|a|=fc 

To state the (analytical and geometrical) conditions under which the problem (|1.2[) . formulated 
as above, is well-posed, we need one final piece of terminology. By the infinitesimal mean oscillation 
of a function F G L\(Q) we shall understand the quantity 

{F} sc(n) ■= limsup ( sup / / F(X) - F{Y) dXdY) , (1.10) 

e^o JB e nn JB E nn J 

where {B e }^i stands for the set of arbitrary balls centered at points of Q and of radius e, and 
the barred integral is the mean value. In a similar fashion, the infinitesimal mean oscillation of a 
function / G Li(<9Q) is defined by 

{f}osc(dn) ■= hmsup sup f f f(X)-f(Y) da x da Y ), (1.11) 



e^o \{B E } gn JB s non JB e ndU 
where {B e }gfi is the collection of n-dimensional balls with centers on d£l and of radius e. 

Our main result reads as follows; see also Theorem Id 101 for a more general version. 

Theorem 1.1 In the above setting, for each p G (1, oo) and s G (0,1), the problem (|1.2|) with 
boundary data as in M.9\) has a unique solution U for which 11. fy) holds with a = 1 — s — 1/p 
provided the coefficient matrices A a p and the exterior normal vector u to dQ satisfy 

Mosc(an)+ Yl {4tf}osc(n) < Cs(l-s){pp' + s- 1 (l-s)- 1 ^ (1.12) 

\a\ = \f3\=m 

where p' = p/{p — 1) is the conjugate exponent of p. Above, C is a sufficiently small constant which 
depends on kq, k% and the Lipschitz constant of£l, and is independent of p and s. Furthermore, 
the bound (jl.121) can be improved for second order operators, i.e. when m = 1, when the factor 
s(l — s) in (|1.12|) can be removed. 



Let BMO and VMO stand, respectively, for the John-Nirenberg space of functions of bounded 
mean oscillations and the Sarason space of functions of vanishing mean oscillations (considered 
either on f2 or on d£l). Since for an arbitrary function F we have (with the dependence on the 
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domain dropped) {-Fjosc < 2 dist (F, VMO) where the distance is taken in BMO, the smallness 
condition (|1.12|) in Theorem II, II is satisfied if 



dist (y, VMO {dty) + ^ dist {A a p, VMO {(!))< C s{l - s)(pp' + S -\l - s)- 1 ^ . (1.13) 

|a|=|j9|=m 

In particular, this is trivially the case when v G VMO(df2) and the A^/j's belong to VMO(f2), 
irrespective of p, s, Kq, K\ and the Lipschitz constant of Q. 

While the Lipschitz character of a domain 0, controls the infinitesimal mean oscillation of its 
unit normal, the inequality in the opposite direction is false in general, as seen by considering 
O := {(x,y) 6R 2 : y > tp £ (x)} with ip e (x) := x sin (e log Indeed, a simple calculation gives 

IKIIbmo(ir) < Ce, yet \\<p' £ \\ Loc (R) ~ 1 uniformly for e G (0,1/2). 

An essentially equivalent reformulation of (|1.2j) is 

D a (A a/3 (X)D /3 U) = in 17, 

\a\ = \/3\=m (1-14) 

Tr [DT W] = 5 7 on «9Q, < m - 1, 

where U satisfies ()1.6|l and 

5 := {2 7 }| 7 |< m -i G i? p m ~ 1+s (^), (1.15) 

though an advantage of the classical formulation (|1.2[) is that the number of the data is minimal. 
For a domain Q C M 2 of class C r , r > ^, and for constant coefficient operators, the Dirichlet 
problem (|1.14|) has been considered by S. Agmon in [2] where he proved that there exists a unique 
solution U G C m ~ 1+S (f2), < s < r, whenever g~ = D^V\qq , I7I < m — 1, for some function 
V G C m " 1+S (^). See also [3j for a related version. 

The innovation that allows us to consider, for the first time, boundary data in Besov spaces as 
in (|1.9|) and (|1.15j) . is the systematic use of weighted Sobolev spaces such as those associated with 
the norm in (|1.6|) . In relation to the standard Besov scale in M n , we would like to point out that, 
thanks to Theorem 4.1 in [25] on the one hand, and Theorem 1.4.2.4 and Theorem 1.4.4.4 in |20j 
on the other, we have 

a = 1 - s - \ g (0, 1 - i/p) w^ a (n) B™~ 1+S+1/P (n), 

(1.16) 

a = 1 - s - J G (-l/p,0) i?; n - 1+s+1/p (0) ^ p m ' a (0). 

Of course, W p ' (O) is just a classical Sobolev space when a = 0. 

Remarkably, the classical trace theory for unweighted Sobolev spaces turns out to have a most 
satisfactory analogue in this weighted context; for the upper half-space this has been worked out 
by S.V. Uspenskh in [49J, a paper preceded by the significant work of E. Gagliardo in |18j in the 
unweighted case. As a consequence, we note that under the assumptions made in Theorem ll.il 



\\Tr[D a U]\\ B s m ~ I [ PiXy^' 1 \D a U(XW dx\ , (1.17) 

|a[<m-l \0<\a\<m Q 
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uniformly in U satisfying C(X, Dx)U = in S7. The estimate (|1.17|) can be viewed as a far-reaching 
generalization of a well-known characterization of the membership of a function to a Besov space 
in W'- 1 in terms of weighted Sobolev norm estimates for its harmonic extension to M" (see, e.g., 
Proposition 7' on p. 151 of |47| ) . 

Theorem 1.1 is new even in the case when m = 1 and A a p € C' x ' (i.e., for second order, constant 
coefficient systems) and provides a complete answer to the issue of well-posedness of the problem 
(jl.2j) . (|l,9j) . (|1.6|) in the sense that the small mean oscillation condition, depending on p and s, is 
in the nature of best possible if one insists of allowing arbitrary indices p and s in (|1,9|) . This can 
be seen by considering the following Dirichlet problem for the Laplacian in a domain f2 C R n : 

' AU = in n, 

< TiU = g£ B s p (dn), (1.18) 

k D a U G L p (n, p(X)P^ 1 -^- 1 dX), Va : \a\ < 1. 

It has long been known that, already in the case when dQ exhibits one cone- like singularity, the 
well-posedness of (|1.18|) prevents the indices (s, 1/p) from taking arbitrary values in (0, 1) x (0, 1). 
At a more sophisticated level, the work of D. Jerison and C.Kenig in [25J shows that (|1.18j) is 
well-posed in an arbitrary, given Lipschitz domain Q if and only if the point (s, 1/p) belongs to a 
certain open region IZq C (0, 1) x (0, 1), determined exclusively by the geometry of the domain £1 
(cf. [25j for more details). Let us also mention here that, even when d£l is smooth and m = I = 1, 
a well-known example due to N.Meyers (cf. |35j) shows that the well-posedness of (|1,2[) in the 
class of operators with bounded, measurable coefficients confines p to a small neighborhood of 2. 

Broadly speaking, there are two types of questions pertaining to the well-posedness of the 
Dirichlet problem in a Lipschitz domain 0, for a divergence form, elliptic system (|1 . If) of order 2m 
with boundary data in Besov spaces. 

Question I. Granted that the coefficients of C exhibit a certain amount of smoothness, identifying 
the Besov spaces for which this boundary value problem is well-posed. 

Question II. Alternatively, for a given Besov space characterize the class of Lipschitz domains Q 
and elliptic operators C for which the aforementioned boundary value problem is well-posed. 

These, as well as other related issues, have been a driving force behind many exciting, recent 
developments in partial differential equations and allied fields. Ample evidence of their impact 
can be found in C. Kenig's excellent account [27\ which describes the state of the art in this field 
of research up to mid 1990's, with a particular emphasis on the role played by harmonic analysis 
techniques. 

One generic problem which falls under the scope of Question I is to determine the optimal 
scale of spaces on which the Dirichlet problem for an elliptic system of order 2m is solvable in an 
arbitrary Lipschitz domain Q in M. n . The most basic case, that of the constant coefficient Laplacian 
in arbitrary Lipschitz domains in W 1 , is now well-understood thanks to the work of B. Dahlberg and 
C. Kenig [14], in the case of L p -data, and D. Jerison and C. Kenig |25j, in the case of Besov data. 
The case of (|1.18|) for boundary data exhibiting higher regularity (i.e., s > 1) has been recently 
dealt with by V. Maz'ya and T. Shaposhnikova in [34J where nearly optimal smoothness conditions 
for d£l are found in terms of the properties of v as a Sobolev space multiplier. Generalizations of 
(TTT81 to the case of variable-coefficient, single, second order elliptic equations have been obtained 
by M. Mitrea and M. Taylor in [31], GUI, |30|. 
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In spite of substantial progress in recent years, there remain many basic open questions, partic- 
ularly for / > 1 and/or m > 1 (corresponding to genuine systems and/or higher order equations), 
even in the case of constant coefficient operators in Lipschitz domains. In this context, one sig- 
nificant problem (as mentioned in, e.g., [16J) is to determine the sharp range of p's for which 
the Dirichlet problem for elliptic systems with L p -boundary data is well-posed. In [42] . J.Pipher 
and G. Verchota have developed a L p -theory for real, constant coefficient, higher order systems 
L = Yl\ a \=2m A a D a when p is near 2, i.e. 2 — e < p < 2 + e with e > depending on the Lipschitz 
character of f2 but this range is not optimal. Recently, more progress for the biharmonic equation 
and for general constant coefficient, second order systems with real coefficients, which are elliptic in 
the sense of Legendre-Hadamard was made by Z. Shen in [44], where he further extended the range 
of p's from (2 — e, 2 + e) to (2 — e, 2 ^~ 3 ^ +e) for a general Lipschitz domain Q in K n , n > 4, where as 
before e = e{dVL) > 0. Let us also mention here the work |T] of V. Adolfsson and J. Pipher who have 
dealt with the Dirichlet problem for the biharmonic operator in arbitrary Lipschitz domains and 
with data in Besov spaces, |51 | where G. Verchota formulates and solves a Neumann-type problem 
for the bi-Laplacian in Lipschitz domains and with boundary data in L2, |37| where the authors 
treat the Dirichlet problem for variable coefficient symmetric, real, elliptic systems of second order 
in an arbitrary Lipschitz domain £1 and with boundary data in Bp(d£l), when 2 — e < p < 2 + e 
and < s < 1, as well as the paper [28J by V. Kozlov and V. Maz'ya, which contains an explicit 
description of the asymptotic behavior of null-solutions of constant coefficient, higher order, elliptic 
operators near points on the boundary of a domain with a sufficiently small Lipschitz constant. 

A successful strategy for dealing with Question II consists of formulating and solving the ana- 
logue of the original problem in a standard case, typically when $7 = Wf and C has constant 
coefficients, and then deviating from this most standard setting by allowing perturbations of a 
certain magnitude. A paradigm result in this regard, going back to the work of Agmon, Doughs, 
Nirenberg and Solonnikov in the 50's and 60's is that the Dirichlet problem is solvable in the context 
of Sobolev-Besov spaces if d£l is sufficiently smooth and if C has continuous coefficients. The latter 
requirement is an artifact of the method of proof (based on Korn's trick of freezing the coefficients) 
which requires measuring the size of the oscillations of the coefficients in a pointwise sense (as 
opposed to integral sense, as in (jl.lUJl ). For a version of Question II, corresponding to boundary 
data of higher regularity, optimal results have been obtained by V. Maz'ya and T. Shaposhnikova 
in [32 J. In this context, the natural language for describing the smoothness of the domain f2 is that 
of Sobolev space multipliers. 

While the study of boundary value problems in a domain £1 C R n for elliptic differential opera- 
tors with discontinuous coefficients goes a long way back (for instance, C. Miranda has considered in 
[36] operators with coefficients in the Sobolev space W^), a lot of attention has been devoted lately 
to the class of operators with coefficients in VMO (it is worth pointing out here that W\ =— > VMO 
on Lipschitz subdomains of M. n ). Much of the impetus for the recent surge of interest in this 
particular line of work stems from an observation made by F. Chiarenza, M. Frasca and P. Longo 
in the early 1990's. More specifically, while investigating interior estimates for the solution of a 
scalar, second-order elliptic differential equation of the form CIA = F, these authors have noticed 
in |10] that IA can be related to F via a potential theoretic representation formula in which the 
residual terms are commutators between operators of C alder on- Zygmund type, on the one hand, 
and operators of multiplication by the coefficients of C, on the other hand. This made it possible 
to control these terms by invoking the commutator estimate of Coifman-Rochberg- Weiss ([12]). 
Various partial extensions of this result can be found in [5], [5], |7|, |19|. |21|. [43], and 

the references therein. Here we would just like to mention that, in the whole Euclidean space, a 
different approach (based on estimates for the Riesz transforms) has been devised by T. Iwaniec 
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and C. Sbordone in |24j. 

Compared to the aforementioned works, our approach is more akin to that of F. Chiarenza 
and collaborators ([10J, [llj), though there are fundamental differences between solving boundary 
problems for higher order and for second order operators. One difficulty inherently linked with the 
case m > 1 arises from the way the norm in (jl.fi j) behaves under a change of variables x : £1 = 
{(X', X n ) : X n > ip(X')} — > R™ designed to flatten the Lipschitz surface d£l. When m = 1, a simple 
bi-Lipschitz changes of variables such as the inverse of the map R" 3 (X' , X n ) i— > (X' , (p(X')+X n ) € 
O will do, but matters are considerable more subtle in the case m > 1. In this latter situation, 
we employ a special global flattening map first introduced by J. Necas (in a different context; 
cf. p. 188 in [41 J) and then independently rediscovered and/or further adapted to new settings by 
several authors, including V. Maz'ya and T. Shaposhnikova in |32j, B. Dahlberg, C. Kenig J. Pipher, 
E. Stein and G. Verchota (cf. | 13| and the discussion in |15|). and S. Hofmann and J. Lewis in [22\. 
Our main novel contribution in this regard is adapting this circle of ideas to the context when one 
seeks pointwise estimates for higher order derivatives of x and A := x _1 in terms of [V(/9] BMO (r™- 1 )- 

Another ingredient of independent interest is deriving estimates for D^DyG(x, y) where G is 
the Green function associated with a constant (complex) coefficient system L(D) of order 2m in 
the upper half space, which are sufficiently well-suited for deriving commutator estimates in the 
spirit of [12J. The methods employed in earlier work are largely based on explicit representation 
formulas for G(x, y) and, hence, cannot be adapted easily to the case of non-symmetric, complex 
coefficient, higher order systems. By way of contrast, our approach consists of proving directly that 
the residual part R(x,y) := G(x,y) — $(x — y), where $ is a fundamental solution for L(D), has 
the property that D^DyR(x,y) is a Hardy-type kernel whenever |a| = \j3\ = m. 

The layout of the paper is as follows. Section 2 contains estimates for the Green function in the 
upper-half space. Section 3 deals with integral operators (of Calderon-Zygmund and Hardy type) as 
well as commutator estimates on weighted Lebesgue spaces. In the last part of this section we also 
revisit Gagliardo's extension operator and establish estimates in the context of BMO. Section 4 
contains a discussion of the Dirichlet problem for higher order, variable coefficient elliptic systems 
in the upper-half space. Then the adjustments necessary to treat the case of an unbounded domain 
lying above the graph of a Lipschitz function are presented in Section 5. Finally, in Section 6, we 
explain how to handle the case of a bounded Lipschitz domain, and state and prove Theorem 16. 1UI 
(from which Theorem 11.11 follows ) . This section also contains further complements and extensions 
of our main result. 



2 Green's matrix estimates in the half-space 

In this section we prove a key estimate for derivatives of Green's matrix associated with the Dirichlet 
problem for homogeneous, higher-order constant coefficient elliptic systems in the half-space R™ . 

2.1 Statement of the main result 

Let L{D X ) be a matrix-valued differential operator 




X 1 



(2.1) 



a\=2m 
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where the A a 's are constant I x I matrices with complex entries. Throughout the paper, D x := 
^-\a\QaiQa2 . . . Qa„ j£ Q _ ( aij qj 2) q,^) g pj™. Here and elsewhere, N stands for the collection of 
all positive integers and No := N U {0}. 

We assume that L is strongly elliptic, i.e. there exists k > such that Y^\ a \=m ll^-allc !x! — K 
and 

»<^(e)»/,»?>C'>«le| 2m Wc'. V£gR", VrjeC. (2.2) 
In what follows, in order to simplify notations, we shall denote the norms in different finite- 
dimensional real Euclidean spaces by | ■ | irrespective of their dimensions. Also, quite frequently, 
we shall make no notational distinction between a space of scalar functions, call it X, and the space 
of vector-valued functions (of a fixed, finite dimension) whose components are in X. 

We denote by F(x) a fundamental matrix of the operator L(D X ), i.e. a I x I matrix solution of 
the system 

L{D x )F(x) = 5{x)k inR n , (2.3) 

where l\ is the I x I identity matrix and 5 is the Dirac function. 
Consider the Dirichlet problem 

f L(D x )u = f inR", 

{ (2-4) 
{ Tr [&>u/dxi] = fj j = 0, 1, . . . , to - 1, on R™ -1 

where R™ := {x = (x',x n ) : x' 6 R™" 1 , x n > 0} and Tr is the boundary trace operator. Hereafter, 
we shall identify cM™ with R n_1 in a canonical fashion. 

For each y' G R n_1 we introduce the Poisson matrices Pq, . . . , P m -i for problem (|2,4|) . i.e. the 
solutions of the boundary- value problems 

[ L(D x )P J (x,y') = 0I l onR™, 

g^Pj) ( 0), y' ) = S jk 5(x' - y% for x' G R n_1 , < fc < m - 1, 

where djfc is the usual Kronecker symbol and < j < m — 1. The matrix- valued function Pj(x,0') 
is positive homogeneous of degree j + 1 — n, i.e. 

Pj(x,0') = \x\ j+1 ~ n Pj(x/\x\,0'), xGR n , (2.6) 

where 0' denotes the origin of R n_1 . The restriction of Pj(;Q') to the upper half-sphere 5? is 
smooth and vanishes on the equator along with all of its derivatives up to order m — 1 (see for 
example, §10.3 in [30j). Hence, 

m 

\\Pi{x,tf)\\ & *i<C (2.7) 

and, consequently, 

\\Pi{xMc>«<C \ x _y*»^ x ^ XGR", y' G R n_1 . (2.8) 

By G(x,y) we shall denote the Green's matrix of the problem (|2.4j) . i.e. the unique solution of 
the boundary- value problem 
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{ L(D x )G(x,y) = 6(x-y)I l for x G R n , 
d j \ 

-G ) ((x',0),2/) =07/ for x' el 



ra-1 



where y £ M™ is regarded as a parameter. 
We now introduce the matrix 



< j < m — 1, 



R(x,y) :=F(x-y)-G(x,y), x,y cR^, 
so that, for each fixed y € M™ , 



(2.9) 



(2.10) 



{ L(D x )R(x,y) = 
i?W,0),y) 



for rc € 



9* 



-F ) ((a/, 0) — for x' € 



»n- 1 



< j < m — 1. 



(2.11) 



.dxi J \dxi 
Our goal is to prove the following result. 

Theorem 2.1 For all multi-indices a, (3 of length m 

\\D^R(x,y)\\ c i x i<C\x-yr n , 
for x, y € , where y := (y', — y n ) is the reflection of the point y 6 with respect to 51 



(2.12) 



In the proof of Theorem 12.11 we distinguish two cases, n > 2m and n < 2m, which we shall 
treat separately. Our argument pertaining to the situation when n > 2m is based on a lemma to 
be proved in the subsection below. 

2.2 Estimate for a parameter dependent integral 

As a preamble to the proof of Theorem 12.11 here we dispense with the following technical result. 

Lemma 2.2 Let a and b be two non-negative numbers and assume that £ 6 WL N . Then for every 
e > and < 5 < N there exists a constant c(N, e,5) > such that 



I 



< 



c(N,s,S) 



(|7?| + a) N +z(\ V - C| + b) N ~ s ~ a £ (|C| + a + b) N ~ s ' 



(2.13) 



Proof. Write J = J\ + J2 where J stands for the integral in the left side of ()2.13|) . whereas J\ 
and J2 denote the integrals obtained by splitting the domain of integration in J into B a = {rj G 
R N : \r]\ < a} and W 1 \ B a , respectively. If |C| < 2a, then 



Jl <a 



-N-e 



d,T] 



< ca 



-N-e 



Hence 



B a (\V~(\ + b) N - s 

-N-e „N /uN-S 



B 4a 



N-S - 



Jl < 



ca - -a'yb 1 
ca- N - £+6 



if a < b, 
if a > b, 



(2.14) 



(2.15) 
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so that, in particular, 

|C| < 2a => Ji < ca- £ (\C\ + a + 6)^. (2.16) 
Let us now assume that |C| > 2a. Then 

Jl £ L s + " + 6)S " (2 - 17> 

which is of the right order. As for J2, we write 



j2 - L\ Ba \v\ N n\v-\\+b) N - s ~ J%1 + j2 > 2 - (2 ' 18) 

where ^2,1, ^2,2 are obtained by splitting the domain of integration in the above integral into the 
set {q : \q\ > max{a, 2|£|}} and its complement in R n \B a . We have 

f drj f drj 

A ~ 4|>max K6 ,2|C|} \v\ N+e M+b) N ~ S + A >k; |> max{aA 2|C|} \v\ N+E M + b) N ~ S 

.( [ + [ 



< C 



drj 

7 |>max{a,6,2|C|} 



(a + b+\(\) N+£ - 5 at(a + b+\(\) 



N-S 



~ at(a + b+ \C\) N - g - (2 " 19) 
There remains to estimate the integral 

j2 > 2 = L m \s a \v\ N+£ (\v- 1 (\+b) N - 5 = J $ + (2 - 20) 

where J 2 ^2 an d J22 are obtained by splitting the domain of integration in J 2 ^ 2 into B\^u 2 \B a and 
its complement (relative to B 2 \Q\\B a ). On the one hand, 



^2 2 < 

' ~ US I -r ") ' JB M/2 \B a l'/r ' " 

r(2) 



On the other hand, whenever \(\ > a/2, the integral J\ \ 2 ', which extends over all 77's such that 
|??| > a, 2\(\ > \r]\ > |C|/2, can be estimated as 

J2 ' 2 - icr +£ J Bm \B a (\v-c\+ b) N - s - icr +£ J Bm m + 
< ^—(l d A + r « 



icr +£ \7i^Ki m+b) N - 5 yiei<4ici +&) Ar -V 



Consequently, 
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(2) cmindCl.ft}^ , . 

^2,2 ^ ^ N+£bN S ■ [Z.ZZ) 

Using [£| > a/2 and the obvious inequality 

minllCl^}^ max{\(\,b N - 5 } < \(\ N b N ~ 5 (2.23) 

we arrive at 

< ca- £ (|C| + a + 6)^. (2.24) 

The estimate (|2.24j) . along with 1)2. 21j) and 1)2.20)) . gives the upper bound ca _e (|£| + a + fe) 5_Ar for 
1/2,2- Combining this with (|2.19|) we obtain the same majorant for J2 which, together with a similar 
result for J\ already obtained leads to 1)2.13)1 . The proof of the lemma is therefore complete. □ 

2.3 Proof of Theorem 12.11 for n > 2m 

In the case when n > 2m there exists a unique fundamental matrix F(x) for the the operator (|2,lj) 
which is positive homogeneous of degree 2m — n. We shall use the integral representation formula 

R(x,y) = R (x,y) + ... + R m „ 1 (x,y), y E M+, (2.25) 

where R(x,y) has been introduced in (|2.1UI) and, with Pj as in Q2.5)l . we set 



Rj(x,y) :-- 



[ Pj(x,£') (^jF) ((g,0)-y)dZ', 0<j<m-l. (2.26) 

JM"- 1 \OXn / 



Then, thanks to (|2.7|) we have 



|# 7 -(x,y)|Ux* < C I -. 7—^7— r^- , , 5 p . (2.27) 



Next, putting 



iV = n — 1 , a = cc n , 
s = m-j , 6 = y„, 
(5 = 2m — j — 1 , £ = y' — x', 

in the formulation of Lemma 12.21 we obtain from 1)2.27)1 

\\Rj(x,y)\\ cM < -fi-, nA _ CX l_ x n _ 2m+7 - , 0<j<m-l. (2.28) 

(\y' - x'\ + x n + y n ) n 2m+ i 



Summing up over j = 0, . . . , m — 1 gives, by virtue of (|2.25j) . the estimate 

\\R(x,y)\\ clxl <C\x-y\ 2m ~ n , x,yeR n + . (2.29) 

In order to obtain pointwise estimates for derivatives of R(x,y), we make use of the following 
local estimate for a solution of problem (|2.4|) with / = 0. Recall that W£ stands for the classical 
Lp-based Sobolev space of order s. The particle loc is used to brand the local versions of these (and 
other similar) spaces. 
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Lemma 2.3 [see [1]] Let ( and £o be functions in Co°(R n ) such that Qq = 1 in a neighbor- 
hood of suppC,. Then the solution u G W|"(R™ , loc) of problem \2.J$ with / = and ipj G 

Wp +1 3 ' 1 ' p (R n ~ 1 , loc), where k >m and p G (l,oo) ; belongs to W k+1 (W\,loc) and satisfies the 
estimate 

TO— 1 

IMI^+i^n) < C llCoPill^+i-i-Vi.pn-1) + I|Cou||l„(r»)) 5 (2-30) 
i=o 

where C is independent of u and <pj . 

Let B(x,r) denote the ball of radius r > centered at x. 

Corollary 2.4 Assume that u G VF™(R™, Zoc) is a solution of problem \2.J$ with / = and 
ifj G C fe+1_ -'(M n_1 , Zoc). Then for any z G and p > 



TO— 1 fe+1— j 

sup |V fe n| < C(p- k sup |u| +V V p s+j ~ k sup |V>j|), (2.31) 

R£n.B(.z,p) ^ WlnB(z,2p) _g S=Q M"- 1 nB(2,2p) ' 

where V' s is the gradient of order s in I"" 1 . Here C is a constant independent of p, z, u and fj. 

Proof. Given the dilation invariant nature of the estimate we seek, it suffices to assume that p = 1. 
Given <f> G C fe+1 ^'( IR?1 ~ 1 ) supported in W 1 " 1 n B(z, 2), we observe that, for a suitable 9 G (0, 1), 



k+l-j 

<C V sup |V>|. (2.32) 

Also, if p > n, 

sup|V fc u| < C\\v\\ w k+i (Rn) , (2.33) 

by virtue of the classical Sobolev inequality. Combining 1)2. 32|) . 1)2. 33 j) with Lemma 12.31 now readily 
gives (E3U). □ 

Given x, y G R" , set p := \x — y\/5 and pick z G <9R™ such that |sg — z\ = p/2. It follows that for 
any w G W±C\B(z, 2p) we have \x—y\ < \x— z\ + \z— w\ + \w— y\ < p/2+2p+\w— y \ < \x—y\/2+\w—y\. 
Consequently, \x — y\/2 < \w — y\ for every w G R™ D B(z,2p), so that, ultimately, 

sup \\V' v F{w-y)\\ cM <- _^_ 2m+fc , (2.34) 

MjeR"-- 1 nB(«,2p) l x ~~ y| 

for each v G No- Granted 1)2.29)1 and our choice of p, we altogether obtain that 

\\D^R(x,y)\\ cM <C k \x-y\ 2m - n - k , x,y £ R™ , (2.35) 
for each multi- index a G Ng of length k. 
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In the following two formulas, it will be convenient to use the notation Rc for the matrix R 
associated with the operator C{D X ) as in ()2.1U|) . By Green's formula 



Rc(y,x) = R c *(x,y) 



x,y G 



where the superscript star indicates adjunction. 

In order to estimate mixed partial derivatives, we observe that (|2.36f) entails 

(D§R c )(x,y)= [(DPR c *)(y,x)Y 
and remark that C* has properties similar to C. This, in concert with (|2.35j) and 



\x — y\ = \x — y\, 



x,y G 



yields 



|2m— n— 



\\DpR(x,y)\\ cM <C p \x-y\ 
Let us also point out that by formally differentiating ()2.11|) with respect to y we obtain 



(2.36) 



(2.37) 

(2.38) 
(2.39) 



£{D x )[DgR c (x,y)]=0 



for x G R n , 



93 D?r) ((x',0),y) = (^-(-DfF^j ((x',0)-y), x> G 



dx 



j y 



I r- TTjn— 1 



dxl 



< j < m — 1. 



(2.40) 



With (|2,39|) and ()2.4()|) in place of (|2.29|) and ()2.11|) , respectively, we can now run the same program 
as above and obtain the estimate 



i2m-n-\a\-\l3\ 



(2.41) 



\\DZDPR(x,y)\\ cM <C a/3 \x-y 
for all multi-indices a and (5. 

2.4 Proof of Theorem 12.11 for n < 2m 

When n < 2m we shall use the method of descent. To get started, fix an integer iV such that N > 2m 



and let (x,z) i— > Q(x,y,z — £) denote the Green matrix with singularity at (y,C) G 
the Dirichlet problem for the operator C{D X ) + (— A 2 ) m in the iV-dimensional half-space 



pn w Tu>N— n 



of 



:= {(x, z): z G R N ~ n , x = (x' , x n ), x' G R n_1 , x n > 0}. (2.42) 

Also, recall that G(x,y) stands for the Green matrix of the problem (|2.4j) . Our immediate goal is 
to establish the following. 



Lemma 2.5 For all multi-indices a and (3 of order m and for all x and y in M™ 

D«D$G(x', y) = j D^Dfeix, y, -Q d(. 



(2.43) 
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Proof. The strategy is to show that 

D^G(x,y)f p (y)dy 



J 



D%DPg(x,y,-Od(f p (y)dy 



(2.44) 



for each fp £ C^°(]R™), from which (|2.43j) clearly follows. To justify (|2.44|) for a fixed, arbitrary 
fp € Co°(lR™), we let u be the unique vector- valued function satisfying D a u € L 2 (R") for all a 
with | a | = m, and such that 



on M n_1 , < j < m - 1. 



C(D x )u = D^ffs in 
\ fd j u\ . . . 

1 U) (I ' 0H0 

It is well-known that for each 7 6 Ng 

\D^u(x)\ < a, |x| m - n - 171 for Id > 1. 



(2.45) 



(2.46) 



This follows, for instance, from Theorem 6.1.4 |29j combined with Theorem 10.3.2 [30J. Also, as 
a consequence of Green's formula, the solution of the problem (j2.45|) satisfies 



D«u(x) = / D«(-D y fG(x, y) fp(y) dy. 



(2.47) 



We shall now derive yet another integral representation formula for D^u in terms of (derivatives 
of) Q which is similar in spirit to ()2.47j) . To get started, we note that since N > 2m the estimate 
(|2.41|) implies 



\D2DPg(x,y,-0\y*i <c(|s-y| + [C[)- 



■N 



(2.48) 



Let us now fix x € R™, p > and introduce a cut-off function H € C°°(M. N n ) which satisfies 
H(z) = 1 for \z\ < 1 and H(z) = for \z\ > 2. We may then write 



u(x) = / G(x, y, -C) H{Clp)D P fp{y) + (- A ? ) m (H (( / p) u{y)) dy d(, 

JR N 1 J 

which further implies 



(2.49) 



D«{-D y f G(x, y, -C) H(C/p) fp(y) dy d( 



<c^ I WDlg{x,v,-O\\ ^\u{y)Dl{H{C/p))\^ (2-50) 



\~/\=m 

By (|2.46|) and (|2.48j) . the expression in the right-hand side of (|2.50|) does not exceed 



n> / d( [ {\ x -y\ + \C\)- N \yr~ n dy 

' p<\(\<2p il™" 1 



< cp 



N—n—m 



(\y\ + p)- N \y\ m ~ n d y = cp- n . 
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This estimate, in concert with (|2.48j) . allows us to obtain, after making p — ► oo, that 

D«u{x)= I [ D«(-D y fg(x,y,-()d(fp(y)dy. (2.51) 

Now d23H) follows readily from this and (f2~T7|) . □ 

Having disposed of Lemma 12.51 we are ready to present the 

End of Proof of Theorem I2.lt Assume that 2m > n and let TV" be again an integer such that 
N > 2m. Denote by F(x,z) the fundamental solution of the operator C{D X ) + (— A z ) m , which is 
positive homogeneous of degree 2m — N and is sing ular at (0,0) el"x R N ~ n . Then the identity 

D°+PF(x)= [ BS+Prfa-Qdt (2.52) 

JR N -" 

can be established by proceeding as in the proof of Lemma 12.51 Combining 1)2.52(1 with Lemma 12.51 
we arrive at 

D%DPR(x',y)= f D^n(x,y,-0d(, (2.53) 
where 1Z(x, y, z) := G(x, y, z) — F(x — y,z). Consequently, 

\\DZDlK(x,y,-0\\ci*t <C(\x-y\ + \C\)' N . (2.54) 

by (|2.35f) with k = and N in place of n. This estimate, together with (j2.53|) . yields (|2.12|) and 
the proof of Theorem I2.1l is therefore complete. □ 



3 Properties of integral operators in a half-space 

In §3.1 and §3.2 we prove estimates for commutators (and certain commutator-like operators) 
between integral operators in W±. and multiplication operators with functions of bounded mean 
oscillations, in weighted Lebesgue spaces on M". Subsection 3.3 contains BMO and pointwise 
estimates for extension operators from R ra_1 onto . 

Throughout this section, given two Banach spaces E, F, we let £(£', F) stand for the space of 
bounded linear operators from E into F, and abbreviate £(E) := £(E,E). Also, given p G [l,oo], 
an open set O C W 1 and a measurable nonnegative function w on O, we let L p (0, w(x) dx) denote 
the usual Lebesgue space of (classes of) functions which are p-th power integrable with respect to 
the weighted measure w(x) dx on O. Finally, following a well-established tradition, A(r) ~ B(r) will 
mean that each quantity is dominated by a fixed multiple of the other, uniformly in the parameter 
r. 



3.1 Kernels with singularity at dM.™ 

Recall L p (R™ , x^ dx) stands for the weighted Lebesgue space of p-th power integrable functions in 
corresponding to the weight w(x) := 4 P , x = (x',x n ) £ M™. 

Proposition 3.1 Let a G R, 1 < p < oo, and assume that Q is a non-negative measurable function 
on {( = ((', ( n ) G R ra_1 x R : Q n > —1}, which also satisfies 
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/ Q(C',Cn-l)C a " 1/P <<00. (3.1) 
JR™ 

Then the operator 

Qf(x):=x- n [ Q(y^)f(y) dy , x=(x',x n )GR n + , (3.2) 

initially defined on functions f G L p (R™) compact support inM™, can 6e extended by continuity 
to an operator acting from L p (M™, dx) into itself, with the norm satisfying 

\\QW P <M, xVdx)) < I Q(C',(n-l)(n a - 1/P d(. (3.3) 
+ JR™ 

Proof. Introducing the new variable ( := {x~ 1 {y l — x'),x~ y n ) G M™, we may write 

\Qf(x)\< f Q(C,(n-l)\f(x' + XnC,XnCn)\d(, VxGM™. (3.4) 

Then, by Minkowski's inequality, 

\\Qf\\L p (Rl, X 7d X ) < [ Q(C',(n-l)([ X^\f(x' + X n C',X n C n )\ P dx) 1/P dC 
+ JR™ V JR™ 7 

= (/ Q(C\Cn-l)Q a ~ 1/P dC)\\f\\ Lpn ^^, (3.5) 

as desired. □ 
Recall that y := (y', -y n ) if y = (y', y n ) G W 1 ' 1 x M. 

Corollary 3.2 Consider 

i?/(x) := / | /(y) dy, ^ = (x', x n ) G R". (3.6) 

|x y\ 



+ 



Then for each 1 < p < oo and eac/i a G (—1/p, 1 — 1 /p) the operator R is bounded from L p (WL, x n dx) 
mto ifee//. Moreover, 

ii on c(n)p 2 c(n) 

\\R\k(L p (R^dx)) < (p a + l)(p(l_ a )_l) = ^T^)' ( } 

where s = 1 — a — 1/p and c(n) is independent of p and a. 
Proof. The result follows from Proposition 13. II with 

._ log (|Cl + 2) 

2(0 - ( | C | 2 + 1)n/2 > (3-8) 

and from the obvious inequality 2\x — y\ 2 > \x — y\ 2 + x^. O 
Let us note here that Corollary 13.21 immediately yields the following. 
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Corollary 3.3 Consider 

Kf(x) := / -^4- dy, x G R» . (3.9) 

Jr™ \x - y\ n 

Then for each 1 < p < oo and a G (— 1/p, 1 — 1/p) the operator K is bounded from , dx) 

into itself. Moreover, 

c(n)p 2 c(n) 
\\K\\ £{LpmjX «v dx)) < {pa + 1){p{1 _ a) _ 1) = ^^-y (3-10) 

where s = 1 — a — 1/p and c(n) is independent of p and a. 

Recall that the barred integral stands for the mean- value (taken in the integral sense). 

Lemma 3.4 Assume that 1 < p < oo, a £ (— 1/p, 1 — 1/p), and recall the operator K introduced 
in (|3.9[) . Further, consider a non-negative, measurable function w defined on WL and fix a family 
of balls T which form a Whitney covering o/R™. Then the norm of wK as an operator from 
L P (R™, x^f dx) into itself is equivalent to 

sup -fw(y) p dy. (3.11) 
B&r jb 

Furthermore, 

c(n) ( f \ i/p 

\\wK\\ z{Lp{R n tX ap dx)) < _ sup( fw(y) p dy) , (3.12) 

where c(n) is independent of w, p, and a. 

Proof. Fix / > and denote by the Euclidean volume of B. Sobolev's embedding theorem 
allows us to write 

n 

ll^/HLw ^ l^l" 1 E \ B \ jP/n \\VjKf\\ P LpiB) , V5 € T. (3.13) 

i=o 

Hence, 



( \x a n w{x){K f)(x)\ p dx < c{n) sup -fw(y) p dy [ x p n a V x{ p \V jK f\ p dx. (3.14) 

«f Bef Jb Jm.™ o<j<i 

Observing that Xn|Vj Kf\ < c(n) Kf and referring to Corollary 13. 31 we arrive at the required upper 
estimate for the norm of wK. The lower estimate is obvious. □ 

We momentarily pause in order to collect some definitions and set up basic notation pertaining 
to functions with bounded mean oscillations. Let / be a locally integrable function defined on W 1 
and define the seminorm 

[/]bmo(R«) : = sup f f(x) - f f(y) dy dx, (3.15) 

B JB JB 
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where the supremum is taken over all balls B in R n . Next, if / is a locally integrable function 
defined on R™ , we set 

[/]bmo(R") : = SU P f f( x ) - f f(y) d V dx, ( 3 - 16 ) 
{B} JsnR™ JBnR™ 

where, this time, the supremum is taken over the collection {B} of all balls B with centers in R™ . 
Then the following inequalities are straightforward 

[/]bmo(R") < sup f f f(x) - f(y) dxdy < 2 [/j B MO(i")- (3-17) 
{B} J BnR^_ JBrarj. 

We also record here the equivalence relation 

[/]BMO(lf) ~ [ Ext /]BMO(R")i (3.18) 

where Ext / is the extension of / onto R n as an even function in x n . 

Finally, by BMO(R") we denote the collection of equivalence classes, modulo constants, of 
functions / on R" for which [/]bmO(R") < °°- 

Proposition 3.5 Let b G BMO(R") and consider the operator 

\b(x)-b(y)\ 



Tf(x) := [ 



\x — y\ n 

Then for each p G (1, oo) and a G (— 1/p, 1 — 1/p) 



f(y)dy, xeR" (3.19) 



T : L P (R™ , x a P dx) — ► L P (R™ , x%> dx) (3.20) 
is a well-defined, bounded operator with 

\\ T h(L p (Ri,x7dx)) < s nJ s \ [ & ]bmo(r^)» ( 3 - 21 ) 
where c(n) is a constant which depends only on n. 
Proof. Given x G R™ and r > 0, we shall use the abbreviations 

b r (x) := / b(y) dy, D r (x) := \b(x) - b r (x)\, (3.22) 

and make use of the integral operator 

Sf(x) := / Dlx ~^J f(y) dy, x G R™ , (3.23) 

Jr™ \ x — y\ 

as well as its adjoint S*. Clearly, for each nonnegative, measurable function / on R™ and each 

x G R%, 

Tf(x) < Sf(x) + S*f(x) + [ 'W*)-_Wy)l f{y)dy 

jr™ \ x ~ y\ 

< Sf(x) + S*f(x) + c(n)[b] BMO{R n ) Kf(x), (3.24) 
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where K has been introduced in ()3.9j) . Making use of Corollary 13.31 we need to estimate only the 
norm of S. Obviously, 

f \b T (x) — b\~..-.\(x)\ 
Sf(x)<D Xn (x)Kf(x) + L*±{ £j4U}.f(y) d y. (3.25) 

jr™ f ~~ y\ 

Setting r = \x — y\ and p = x n in the standard inequality 

\b p (x) - b r (x)\ < c(n) log (- + l) [6]bmo(K«), ( 3 - 26 ) 
where r > p (cf., e.g., p. 176 in |32j, or p. 206 in [48J ) , we arrive at 

Sf(x) < D Xn (x)Kf(x) + c(n) [b] BMO (R V Rf(x), (3.27) 

where R is defined in (|3.6|) . Let T be a Whitney covering of Wf with open balls. For an arbitrary 
B G T , denote by 5 the radius of B. By Lemma 13.41 with w(x) := D Xn (x), the norm of the operator 
D Xn (x)K does not exceed 



supf / \D Xn (x)\ p dx) P < c(n) sup ( f \b(x) - b s (x) \ p dx) P + c(n) [6]bmo(ki) 
BeF^JB ' b^t^Jb ' 

< c(n) [6]bmo(r^)> ( 3 - 28 ) 

by the John-Nirenberg inequality. Here we have also used the triangle inequality and the estimate 
(|3.26|) in order to replace b Xn (x) in the definition of D Xn (x) by bg(x). The intervening logarithmic 
factor is bounded independently of x since x n is comparable with 5, uniformly for x S B. With 
this estimate in hand, a reference to Corollary 13.21 gives that 



S : L p (R™ , < p dx) -» L P (IR^, < p dx) boundedly (3.29) 

for each p G (1, oo) and each a G (— 1/p, 1 — 1/p). 

The corresponding estimate for the norm S results as well. By duality, it follows that S* enjoys 
the same property and, hence, the operator T is bounded on L P (W^_, x^f dx) for each p G (l,oo) 
and o G (—1/p, 1 — 1/p), thanks to (J3.24|) and Corollary 13.31 The fact that the operator norm of 
T admits the desired estimate is implicit in the above reasoning and this finishes the proof of the 
proposition. □ 

3.2 Singular integral operators 

We need the analogue of Proposition 13.51 for the class of Mikhlin-Calderon-Zygmund singular inte- 
gral operators. Recall that 

Sf(x)=p.v. [ k(x,x-y)f(y)dy, x G R n , (3.30) 

(where p. v. indicates that the integral is taken in the principal value sense, which means excluding 
balls centered at the singularity and then passing to the limit as the radii shrink to zero), is called 
a Mikhlin-Calderon-Zygmund operator (with a variable coefficient kernel) provided the function 
k : R n x (M. n \ {0}) -» R satisfies: 
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(i) k(x, •) G C°°(R n \ {0}) and, for almost each x G R n , 




max \\D^k(x, z)\\ Loc(R n xS n-i ) < oo, 



(3.31) 



where is the unit sphere in R n ; 

(ii) /c(:r, Az) = \~ n k(x, z) for each z G R n and each A G R, A > 0; 

(iii) J 5 „-i k(x,uj) doj = 0, where da; indicates integration with respect to uj G S n . 

It is well-known that the Mikhlin-Calderon-Zygmund operator 5 and its commutator [S, b] with 
the operator of multiplication by a function b G BMO(R") are bounded operators in L P (R") for 
each 1 < p < oo. The norms of these operators admit the estimates 



where c(n) depends only on n and the quantity in 1)3.31 Jl . The first estimate in 1)3.32)1 goes back 
to the work of A. C alder on and A. Zygmund (cf., e.g., |5], see also the comment on p. 22 of 
[47J regarding the dependence on the parameter p of the constants involved) . The second estimate 
in (|3.32|) was originally proved for convolution type operators by R. Coifman, R. Rochberg and 
G.Weiss in [12J and a standard expansion in spherical harmonics allows to extend this result to 
the case of operators with variable-kernels of the type considered above. 

We are interested in extending (|3.32|) to the weighted case, i.e. when the measure dx on R" is 
replaced by Xn P dx, where 1 < p < oo and a G (—1/p, 1 — l/p). Parenthetically, we wish to point 
out that a G (—l/p, 1 — l/p) corresponds precisely to the range of a's for which w(x) := Xjx* is a 
weight in Muckenhoupt's A p class, and that while in principle this observation can help with the 
goal just stated, we prefer to give a direct, self-contained proof. 

Proposition 3.6 Retain the above conventions and hypotheses. Then the operator S and its com- 
mutator [S,b] with a function b G BMO(R") are bounded when acting from L p (R!f, x^f dx) into 
itself for each p G (1, oo) and a G ( — l/p, 1 — l/p). The norms of these operators satisfy 





(3.32) 




(3.33) 




(3.34) 



Proof. Let Xj be the characteristic function of the layer 2^ 2 < x n < 2 1+ ^ 2 , j = 0, ±1 
Yljez Xj = 2. We then write S as the sum Si + S2, where 



so that 




(3.35) 



\j-k\<3 



The following chain of inequalities is evident 
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\SifW 



Lp(R1 t Xn P dx) 



< 



E / 5 ( E **/ 



fc-j|<3 



j Jk + \k-j\<3 
In concert with the first estimate in ()3.32j) . this entails 



S( E X k 2^ 2 f)( 



p \ i/p 



p \i/p 
dx 



(3.36) 



l^i/II^CRn,^^) < c{n) VP '{Vf (J2 Xk2 ja/2 \f\Y dx 



p \ 1/p 



J + |*-j|<3 

< c(n)pp'(f \f(x)\Pa%>dx) 1/P , 



which is further equivalent to 



IISiII^pCR",^^)) <c(n)pp'. 
Applying the same argument to [5, 6] and referring to ([3.32ft . we arrive at 



(3.37) 
(3.38) 

we arrive at 

II [Si,b] 1 1 £(i/j, (r» , Xn p dx)) < c(n)pp' [6]bmo(R™)- (3-39) 

It remains to obtain the analogues of (|3.38|) and ()3.39j) with 52 in place of Si. One can check 
directly that the modulus of the kernel of 52 does not exceed c(n) \x — y\~ n and that the modulus 
of the kernel of [52,6] is majorized by c(n) |6(x) — 6(y)| \x — y\~ n . Then the desired conclusions 
follow from Corollary 13.31 and Proposition 13.51 □ 

3.3 The Gagliardo extension operator 

Here we shall revisit a certain operator T, extending functions defined on M"" 1 into functions 
defined on M™, first introduced by Gagliardo in |18|. Fix a smooth, radial, decreasing, even, 
non-negative function £ in M n_1 such that £(t) = for \t\ > 1 and 



J ((t)dt = i. 



(A standard choice is £(i) := cexp (—1/(1 — |t| 2 )+) for a suitable c.) Following 



(T(p)(x',x n )= J C(t)tp(x? + x n t) dt, (x',x„)eR™, 



(3.40) 

we then define 
(3.41) 



acting on functions ip from Li(R n ,loc). To get started, we note that 



V x ,(T<p)(x',x n ) 




dx n 



(Ttp)(x',x n ) 



C(t)\7ip(x' + tx n )dt, 
C(t)tVcp(x' + tx n )dt, 



(3.42) 
(3.43) 
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and, hence, we have the estimate 



(^)I!l 00 (R™) < c||Va;/ V9|| Loo ( M n-i). 
Refinements of (|3.44j) are contained in the Lemmas 13.713.81 below. 

Lemma 3.7 (i) For all x 6 and for all multi-indices a with \a\ > 1, 



BMO(R n_1 ) • 



(ii) For all x = (x',x n ) € R%, 

(T(f)( X ) - ip{x') < CX n [V^] BMO (Rn-l ) . 

Proof. Rewriting 1)3.43(1 as 



(3.44) 



(3.45) 



(3.46) 



_d_ 

we obtain 



■(T<p)(x',x n ) = x n 



l-n 



£-x'\£-x' 



V<p(€)-+ Vif{z)dz)di 

x ri. v J\z— x'\<x n 



d 

OX n 

for every non-zero multi- index 7. Furthermore, for i = 1, . . . n — 1, by (|3.42|) 



d 

—V x i(Tip){x) = x 



l-n 
n 



(3.47) 



(3.48) 



/ ^c(— )(V^)-/ V<p(z)dz)d£, (3.49) 

J V x n ' \ J\z-x'\<x n ' 



where di is the differentiation with respect to the i-th component of the argument. Hence once 
again 



d 

D l-Q^T^x'{Ttp){x) < c^ l7l " 1 [V99] BMO ( R n-i), 



(3.50) 



and the estimate claimed in (i) follows. 

Finally, (ii) is a simple consequence of (i) and the fact that (TV)| K «-i = (p. □ 

Remark. In concert with Theorem 2 on p. 62-63 in |47j , formula Q3.42JI yields the pointwise 
estimate 

|V(2V)(jc)| < cM(V<p)(x'), x = (x' , x n ) E R" , (3.51) 

where M. is the classical Hardy-Littlewood maximal function (cf., e.g., Chapter I in |47j). As for 
higher order derivatives, an inspection of the above proof reveals that 



D«(T<p)(x) <cx*-W(V<p)*(x'), (x',x n )eR n , 



(3.52) 



holds for each multi-index a with \a\ > 1, where is the Fefferman-Stein sharp maximal function 

(cf. mi). 
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Lemma 3.8 IfV x ,Lp G BMO^"" 1 ) then V(7» G BMO(R^) and 

[V(TV)] BMO ( R ») < c[V x '^] BMO (K»-i). (3.53) 

Proof. Since (Ttp)(x', x n ) is even with respect to x n , it suffices to estimate [V x (Tip)]BMO(R n )- Let 
Q r denote a cube with side-length r centered at the point 77 = (jj',r) n ) G M n_1 x R. Also let Q' r be 
the projection of Q r on R n_1 . Clearly, 



V x/ (T<p)(x',x n )-V x ,<p(x')=x 1 - n / C 



Suppose that |r/ n | < 2r and write 



(V^O-V^(x'))^. (3.54) 



/ 

Jq. 



V x '(Tcp)(x',x n ) - V x >ip(x') 



dx < cr 2 ~ n I I \V(p{0 -V<p(z)\dzd£. 
JQ' 4r JQ' 4r 

< cr n [V^BMO(R«-i)- (3-55) 



Therefore, for \r] n \ < 2r 



t t \^x'Tip(x) — V y 'Tip(y)\dxdy < 2 + \ V x 'Tip(x) — ~S/(p(x')\ dx 

JQrJQr ■'<,>, 

+ / / \V<p(x')-V<p(y')\dx'dy' 

J Q'r J Q'r 

< c[V(/?] BMO (Rn-i). (3.56) 

Next, consider the case when \rj n \ > 2r and let x and y be arbitrary points in Q r ( r ?)- Then, 

using the generic abbreviation Je '■= + f, we may write 

Je 



\V x 'T(p(x) - V y >Tip(y)\ < 



v^(0 - 



2\Vn\ 



< 



c r 



2|iJnl 



*2|»Jnl 

< c[V<^] BMO (Rn-l). 



Consequently, for \r) n \ > 2r, 

i i \V X /T(p(x) - V y iT(p(y)\dxdy < c [V^ BMO(r -i) 

J Qr Qr 

which, together with (|3.56|) . gives 

[Vx-r^BMOfl") < c[V99] BMO (]In-i). 

This inequality and (|3.48j) . where I7I = 0, imply (|3.53|) . 



(3.57) 



(3.58) 



(3.59) 
□ 
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4 The Dirichlet problem in R" for variable coefficient systems 

4.1 Preliminaries 

For 

1 < p < oo, — < a < 1 and m € N, (4.1) 

we let V^ m,a (IR!L) denote the weighted Sobolev space associated with the norm 



MIVp m ' a (R") : = 



£ / |xlfl-^u(x)r< a da ; ) 1/P (4.2) 



0<|/3|<m 



It is easily proved that C^°(M™ ) is dense in V™' a (W^_). Moreover, by the one-dimensional Hardy's 
inequality (see, for instance, [31], formula (1.3/1)), we have 

IHIv$»-^)^ cs " 1 (E / \D P u{x)\v X i a dx) l/P for«£C„T + ). (4.3) 

\P\=m Ju 1 

The dual of V^'' a (RD will be denoted by V™' a (Wl), where 1/p+l/p' = 1. 
Consider now the operator 

L(x,D x )u:= £ ^(^(xJxW+^-^ti) (4.4) 

0<H,|/3|<m 

where -A Qi g are C Zx '-valued functions in L 00 (M™). We shall use the notation L(x,D x ) for the 
principal part of L(x,D x ), i.e. 

L(x,D x )u: = Yl D a x (A aP (x)D p x u). (4.5) 

\a\ = \f}\=m 

4.2 Solvability and regularity result 

Lemma 4.1 Assume that there exists k = const > such that the coercivity condition 

K /„ E (A a p(x)D^u(x), D a u(x)) c i dx > k J2 P^lli 2 (R«), (4-6) 

^+ |a| = |/3|=m \-y\=m 

holds for all « G C °°(R^), and that 

WAxphcc^D < k -1 - (4-7) 

|a|=]/3|=m 

(i) Lei p G (1, oo) and —l/p<a< 1 — 1/p and suppose that 

g ( 1 _ n ^ IIA^IIl^R") + XT [^a/)]BMO(f!f) - ( 4 - 8 ) 

^ ' \<x\ + \p\<2m \ a \ = \f3\=m 

O<\a\,\0\<m 

where 5 satisfies 
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pp H — jr. r ) ^ < c(n, m, k) (4-9) 



a(l - s 

luzf/i a sufficiently small constant c(n, m, n) and s = 1 — a — 1/p. 
Then the operator 

L = L(x, D x ) : V p m ' a (Rl) — > F p - m ' a (R!J.) (4.10) 

is an isomorphism. 

(ii) Let pi £ (1, oo) and —1/pi < a» < 1 — where i = 1,2. Suppose that \4-yj ) holds with 

Pi and Si = 1 — ai — 1/pi in place of p and s. If u £ Vp™' ai (M+) is sttc/i £/iai Lu € V^ m ' ai (R+) H 

^ m ' a2 (Ky,ttenuGC 2 (^)' ' 

Proof. The fact that the operator (|4,1()|) is continuous is obvious. Also, the existence of a bounded 
inverse L~ l for p = 2 and a = follows from l)4.6j) and (|4.8j) - (|4.9|) with p = 2, a = 0, which allow 
us to implement the Lax-Milgram lemma. 

We shall use the notation L y for the operator L(y,D x ), corresponding to (|4.5j) in which the 
coefficients have been frozen at y € K™, and the notation G y for the solution operator for the 
Dirichlet problem for L y in M™ with homogeneous boundary conditions. Next, given u £ V p m ' a (W^ ), 
set / := Lu € F p ~ m ' a (M^) so that 

L(x,D)u = f in R™, 

&> u , (4.11) 
^(x',Q) = Q onl"" 1 , 0< j <m-l. 

Applying the trick used for the first time in |10| , we may write 



u(x) = (G y f){x)-(G y (L-L y )u)(x)-(Gy{L-L)u){x), x £ R™ . (4.12) 

We desire to use (|4.12j) in order to express u in terms of / (cf. (|4.27() - (|4,28|) below) via integral 
operators whose norms we can control. 

First, we claim that whenever I7I = m, the norm of the operator 



V p m > a (R n + ) Dl(G y (L - L y )u)(x) G L p (R n + , y%> dy) (4.13) 

x=y 



does not exceed 



C (pp' + si l_ s ^ ) Yl [4tf]BMO(R$)- ( 4 - 14 ) 

' \ a \=\/3\=m 

Given the hypotheses under which we operate, the expression (|4.14[) is therefore small if 5 is small. 

In what follows, we denote by G y (x, z) the integral kernel of G y and integrate by parts in order to 
move derivatives of the form D" with \a\ = m from (L—L y )u onto G y (x, z) (the absence of boundary 
terms is due to the fact that G y (x, •) satisfies homogeneous Dirichlet boundary conditions). That 
()4.14|) bounds the norm of (|4.13l) can now be seen by combining Theorem 12.11 with 1)3.21(1 and 
Proposition 13.61 

Let 7 and a be multi-indices with I7I = m, \a\ < m and consider the assignment 
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3*^(d2 [ G y (x,z)D 



m— \a\ 



x=y 



After integrating by parts, the action of this operator can be rewritten in the form 

-1 d \"H C 



Dl 



i dz r , 



(-D z ) a G y (x,z) T a (z)dz 



x=y 



(4.15) 



(4.16) 



where 



Tjz) := < 



\&(z), if \a\ = m, 
^ ]_^ m— \ a \ f°° 



(t-Zn)" 1 -!"'- 1 ^,^ dt, if \<\ <.n>. 



(4.17) 



(m - \a\ - 1)! J Zn v " f m -l°l 
Using Theorem 12. II along with Q3.21JI and Proposition 13.61 we may therefore conclude that 



D2 



-1 \ m ~ l a 



i dz n ' 



x=y 



D z ) a G y (x,z)\T a (z) dz } 



On the other hand, Hardy's inequality gives 



< 



C 



(4.18) 



(4.19) 



II 1 a\\L p {W\.,xl p dx) - ^ _ ll^llip^.x^rfx) 

and, hence, the operator (|4.15|) can be extended from Cq°(IR™ ) as a bounded operator in L p (R™ , x^f dx) 
and its norm is majorized by 



C 



pp + 



l-s\^ 1 s (l-s)J- 
Next, given an arbitrary u £ V^ n ' a (W^), we let \l/ = in (|4.15|) with 

V a/3 (z) := zf- m A ap D"u{z), \a\ + \(3\ < 2m, 
and conclude that the norm of the operator 



V p m > a (R n + ) 3u^Dl (G y (L - L)u){x) € L P (W + , y%> dy) 



does not exceed 



x=y 



^ ^2 \\ A <xp\\Loo(Wt 



PP + 



1-sV* s(l - s) 



(4.20) 

(4.21) 
(4.22) 

(4.23) 



|a| + |3l<2m 
0<|a|,|/9|<m 



It is well-known (cf. (1.1.10/6) on p. 22 of [31J) that any u G Vp l,a (R%) can be represented in 
the form 



u = K{D a u}\ 



(4.24) 
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where K is a linear operator with the property that 



D a K : L P (R™ , < p dar) — » L p (R™ , < p da;) (4.25) 
is bounded for every multi- index a with |a| = m. In particular, by (|4.3|) . 

||if{Z) CT 7i} k | =m || v m,a (R „ ) < Cs^^K^nji^i^H^^n^P^). (4.26) 

At this stage, we transform the identity Q4.12JI in the following fashion. First, we express the 
two u's occurring inside the Green operator G y in the left-hand side of (|4.12|) as in (|4.24jl . Second, 
for each multi-index 7 with I7I = m, we apply D 7 to both sides of (|4.12|) and, finally, set x = y. 
The resulting identity reads 

{DM| 7 |= ra + S{D°u} lal=m = Qf (4.27) 

where Q is a bounded operator from V^~ m ' a (W^_) into L p (R", x^f dx) and 5 is a linear operator 
mapping L P (R™, xff dx) into itself. Furthermore, on account of (j4~T3|) - (j4~TH) . (|4.22j) - (|4.23fl and 
625), we can bound II^I^i^r",^^)) by 

C { pp ' + ^ [^]bmo(K!^) + fl _ v S IK/3||i,°°(R ? ))- (4-28) 

V ^ U| = |/3|= m V ^ M + |,3|<2m 

0<|a|,|/3|<m 

Owing to (|4.8|) - (|4.9|) and with the integral representation formula Q4.27JI and the bound Q4.28JI 
in hand, a Neumann series argument and standard functional analysis allow us to simultaneously 
settle the claims (i) and (ii) in the statement of the lemma. □ 

5 The Dirichlet problem in a special Lipschitz domain 

In this section as well as in subsequent ones, we shall work with an unbounded domain of the form 

G = {X = (X',X n ) eK n :I'e R n -\ X n > <p(X')}, (5.1) 
where if : R n_1 — > R is a Lipschitz function. 

5.1 The space BMO(G) 

The space of functions of bounded mean oscillations in G can be introduced in a similar fashion to 
the case G = W\_ . Specifically, a locally integrable function on G belongs to the space BMO(G) if 

[/]bmo(0 == / f(X) - / f(Y) dY dX < 00, (5.2) 

{B} JBnG JBnG 

where the supremum is taken over all balls B centered at points in G. Much as before, 

[/]bmo(G) ~ sup / / fiX) - f(Y) dXdY. (5.3) 

{B} JBnG JBnG 

This implies the equivalence relation 

1/]bmo(g) ~ [/ A]bmo(r™) (5-4) 
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for each bi-Lipschitz diffeomorphism A of R" onto G. As direct consequences of definitions, we also 
have 

[ ]J /j']bmo(G) < c II/I|J!^(g)[/]bmo(G)> where / = (/i,... ,/tv), (5.5) 
l<j<N 

I/ _1 ]bmo(g) < c II/ _1 ||| 00 (g)[/]bmo(g)- (5-6) 
5.2 A bi-Lipschitz map A : Wl — > G and its inverse 

Let ip : M n_1 — > R be the Lipschitz function whose graph is dG and set M := llV^Hi^Rn-i). Next, 
let T be the extension operator defined as in (|3.41l) and, for a fixed, sufficiently large constant 
C > 0, consider the Lipschitz mapping 

A : R" 3 (x',x n ) i— »• (X',X n ) e G (5.7) 

defined by the equalities 

X':=x', X n :=CMx n + (Tip)(x',x n ) (5.8) 
(see [32J, §6.5.1 and an earlier, less accessible, reference |33|). The Jacobi matrix of A is given by 

V= ( V*>(T<P) CM + d(T^)/dx n ) (5 ' 9) 

where / is the identity (n — 1) x (n— l)-matrix. Since \d(Tip)/dx n \ < cM by (jSHSj) ) it follows that 
det A' > (C - c) M > 0. 

Next, thanks to and (j5~Tj) - (|5~%|) we have 

X n - ppr') ~ x n . (5.10) 

Also, based on (|3.53f) we may write 

]bMO(R™ ) - cfV^BMOCR™- 1 ) ( 5 -H) 

and further, by and (|3~I5|) . 

||Z) a A'(x)|| R „xn < c(M) x-I q I[V^]bmo(R-i), Va : |a| > 1. (5.12) 

Next, by closely mimicking the proof of Proposition 2.6 from [34J it is possible to show the 
existence of the inverse Lipschitz mapping k := A -1 : G — > R" . Owing to (|5.4|) . the inequality 
1)5.11(1 implies 

[A' o x] BMO (G) < cfV^BMOCR"- 1 )- (5.13) 
Furthermore, (|5.12|) is equivalent to 

||(£ a A')(x(X))|[ R „xn < c(M,a)(A n - V 9(A / ))- |q| [V V 9] B mo(r-i), (5-14) 
whenever \a\ > 0. Since x' = (A' o x) _1 we obtain from ()5.6() and (|5.13|) 

Mbmo(G) < cIV^bmoCR"- 1 )- ( 5 - 15 ) 
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On the other hand, using x' = (A' ox) 1 and (|5.14|) one can prove by induction on the order of 
differentiation that 

||Z7V(A0|| R „xn < c(M,a)(X n -^(X / ))- |Q| [V^ BMO ( R n-i ) (5.16) 
for all X G G if \a\ > 0. 

5.3 The space V™> a (G) 

Analogously to V™' a (W^_), we define the weighted Sobolev space V™' a (G) naturally associated with 
the norm 

\Mv^°(G)--=( E / \(X n -^X')p- m D^(XW(X n - i p(X')rdXy /P . (5.17) 

0<N<m 

Replacing the function X n — f(X') by either p(X) := dist (X,dG), or by the so-called regularized 
distance function p Teg (X) (defined as on pp. 170-171 of |47| ) . yields equivalent norms on V™' a (G). 
Based on a standard localization argument involving a cut-off function vanishing near dG (for 
example, take rj(p reg /e) where r/ G C^°(M) satisfies rj(t) = for \t\ < 1 and rj(t) = 1 for \t\ > 2) one 
can show that Cq°(G) is dense in V™ ,a (G). 
Next, we observe that for each U G C^(G), 

Cs\\U\\ vr « {G) <(j2 [ \D^(X)\P(X n -^X')rdx) 1/P < \\U\\ Vp ^ {G) (5.18) 

II J G 
|7|=m 

where, as before, s = 1 — a — 1/p. Indeed, for each multi-index 7 with < I7I < m, the one- 
dimensional Hardy's inequality gives 

/ \(X n - <p(X')p- m D~>U(X)\v (X n - ^{X')f a dX 
JG 

<(C/sfJ2 [ \D a U(X)\P(X n -<p(X')y a dX, (5.19) 

11 " G 
\a\=m 

and the first inequality in (j5.18j) follows readily from it. Also, the second inequality in (|5.18j) is a 
trivial consequence of (|5.17|) . 

Going further, we aim to establish that 

ci IMIy^^n) < \\u o x\\ v ™*(g) < c 2 INly"'"^"), (5.20) 

where c\ and C2 do not depend on p and s, whereas x : G — > Wl is the map introduced in §5.2. 
Clearly, it suffices to prove the upper estimate for ||^°>f||y p m ' ti (G) m (|5.20|) . As a preliminary matter, 
we remark that 

D^u(x(X))) = ((x'*(A)0| =D u)(x(A)) 

n 

l<|r|<| 7 | <r i=l j 
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where 

J2<Tij = 7, > 1, Y^(Wij\ - 1) = l7l - \r\- (5.22) 

In turn, ()5.21jl - (|5.22|) and 1)5. 16j) allow us to conclude that 

\DUu(x(X)))\<c Yl ^-^\D T u{x)l (5.23) 
i<|t|<M 

which, in view of (|5.1U|) . yields the desired conclusion. 
Finally, we set 

V- m > a (G) := (vp~ a {G))*. (5.24) 

where, as usual, p' = p/(p — 1). 

5.4 Solvability and regularity result for the Dirichlet problem in the domain G 

Let us consider the differential operator 

CU = £(X,D X )U = Yl D a {% aP {X)D^U), X € G, (5.25) 
H=|/J|=m 

whose matrix-valued coefficients satisfy 

E ii^iiwo ^ ^- ( 5 - 26 ) 

\a\=\P\=m 

This operator generates the sesquilinear form £(•, •) : V™ ,a {G) x V^ l ~ a {G) — > C where p' is the 
conjugate exponent of p, defined by 



C(U,V):= V [ ^(X)D^U(X), D a V(X))dX. (5.27) 

I., \Z! JG 

We assume that the inequality 



\a\ = \P\=m G 



$t£(U,U)>K Y, W WU WUg) (5-28) 

| 7 |=m 



holds for all hi 6 F 2 m,0 (G). 



Lemma 5.1 (i) Lei p € (l,oo) ; — 1/p < a < 1 — 1/p and s := 1 — a — 1/p. Suppose that 

[V^BMOGR"" 1 ) + ^"^BMOfG) < ^ (5-29) 
|a| = |/3|=m 

where 5 satisfies 

{PP + 1(13^)) ^T^l) < C(n,m,K, ||V^|| ioo(K n-i)) (5.30) 
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with a sufficiently small constant C, independent ofp and s. In the case m = 1 the factor 5 / s(l — s) 
in ()5.3U|) can be replaced by 5. 
Then the operator 



C(X, D X ) : V p m > a (G) — V- m > a (G) (5.31) 

is an isomorphism. 

(ii) Let pi £ (1, oo) and —l/p% < a« < 1 — 1/pi, where % = 1,2. Suppose that Q4.9|) ZioWs imi/i pi 
and si = 1 — at — 1/pi in place of p and s. IftiE V p ™' ai (G) and CU € Vp~ m ' ai (G) n V p ~ 2 m,a2 (G) , 
thenU gC(G). 

Proof. We shall extensively use the flattening mapping A and its inverse studied in §5.2. The 
assertions (i) and (ii) will follow directly from Lemma 14. II as soon as we show that the operator L 
defined in R" by 

L(U o A) := (CU) o X (5.32) 

satisfies all the hypotheses in that lemma. The sesquilinear form corresponding to the operator L 
will be denoted by L(u,v). 

Set u(x) := U(X(x)), v(x) := V(X(x)) and note that the identity (|5.21l) implies 

D0U(X) = ({*»(\(x))t)$ =D u)(x)+ Yl Kp T {x)x\^\D T u{x), (5.33) 

1<M<|/3| 

D a V{X) = {{x'*{X{x))Ot=Dv){x)+ Y K aT (x)x^D T v(x), (5.34) 

i<M<H 

where, thanks to (|5.16|) . the coefficients K^ T satisfy 



ll-^TrlUooOR^.) < cIV^BMOCIR™- 1 )- (5.35) 

Plugging (|5.33f) and (|5.34|) into the definition of C(U, V), we arrive at 

C(U,V) = L (u,v)+ Yl [ (A a /3(x)xt Wl ~ 2m D^u(x), D a v(x))dx, (5.36) 

l<|a|,|fl|<m 



C|c|,|/3|<m •"?>■+ 
\a\ + \/3\<2m 



where 



L (u,v)= V / ((Vi a poX)((x'*oX)OL D u, (^*oX)^ =D v)detX'dx. (5.37) 

\a\ = \/3\=m JM + 

It follows from (|5.33|) - (|5.35j) that the coefficient matrices A a/ 3 obey 

Y, 1 1 A*y8 1 1 < CK ~ l [V^BMOCIR"- 1 )) (5.38) 

l<|a|,|/3|<m 
|a| + |/9|<2m 

where c depends on m, n, and ||Vc/?|| ioo ( R n-i). We can write the form Lq(u,v) as 

Y [ {A aP {x)D p u{x), D a v{x))dx (5.39) 

\ a \ = \f3\= m jR l 
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where the coefficient matrices A a p are given by 

A aP = det\' Y, ^(^'oA)(2L TT oA) 5 (5.40) 

l7l = l T l =m 

for some scalar homogeneous polynomials of the elements of the matrix V(A(x)) with deg P7g = 
2m. In view of (f53 J) - (|5.15[) . 

Y K/?]bmo(R^) < c ( K ~ 1 [ v '^] BMO ( r -i) + Y Mbmo(G)). ( 5 -41) 

M = |/3|=«i |a| = |/3|=m 

where c depends 

By oa 



and, therefore, 



- L (u,u)\ < c5||-u||^ m , 0(R „ ) (5.42) 



KL (n,n) > RC(U,U) - cS\\u\\ 2 vr , 0m) . (5.43) 



Using (|5.28f) and the equivalence 

||W|| y m,o (G) ~ IMIy™' (R+) (5-44) 

(cf. the discussion in §5.3), we arrive at (|4.6|) . Thus, all conditions of Lemma 14.11 hold and the 
result follows. The improvement of 1)5.30(1 for m = 1 mentioned in the statement (i) holds because 
in this case L = Lq. □ 



6 Dirichlet problem in a bounded Lipschitz domain 
6.1 Preliminaries 

Let f2 be a bounded Lipschitz domain in M. n which means (cf. [47J, p. 189) that there exists a finite 
open covering {0j}i<j<iv of <9$7 with the property that, for every j € {1, ...,N}, Oj n 0, coincides 
with the portion of Oj lying in the over-graph of a Lipschitz function tpj : W 1 ^ 1 — > R (where 
M n_1 x R is a new system of coordinates obtained from the original one via a rigid motion). We 
then define the Lipschitz constant of a bounded Lipschitz domain f2 C M n as 

inf (max{||V^ 11^(^-1) : 1 < j < N}) (6.1) 

where the infimum is taken over all possible families {<fj}i<j<N as above. 

It is a classical result that the surface measure da is well-defined and that there exists an 
outward pointing normal vector v at almost every point on d£l. 

We denote by p{X) the distance from X £ W 1 to d£l and, for p, a and m as in (|4.1|1 . introduce 
the weighted Sobolev space V^ m ' a (f2) naturally associated with the norm 

\\u\\v™^) ■■= ( E f \pm^i^{X)rp{xrdxf ,p . m 

Q<\P\<m 

One can check the equivalence of the norms 
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\V p m ' a (Q) ~ HP?eg W lly». (n)' ( 6 - 3 ) 

where p reg (X) stands for the regularized distance from X to dVt (in the sense of Theorem 2, p. 171 
in |47| ) . It is also easily proved that Cg°(fi) is dense in Vp n ' a (Q) and that 



(6.4) 



(2 / l^(x)| p p(x) pa dx) 1/p 

|/3|=m 

uniformly for U G C£°(fi). As in (jS^U), we set 

y- m ' a (^) : =(y p ?'- a (o))*. (6.5) 

Let us fix a Cartesian coordinates system and consider the differential operator 

AU = A{X,D X )U := Yl D a (A aP {X)D /3 U), lefi, (6.6) 

\a\ = \/3\=m 

with measurable Ixl matrix- valued coefficients. The corresponding sesquilinear form will be denoted 
by A(U, V). Similarly to H5.26j) and (|5.28f) we impose the conditions 

IIAtflUcPl) < « -1 (6.7) 

\ct\=\/3\=m 

and 

UA{U,U) >^Y1 U Wl 2 (n) for a11 u e V 2 m '°{G). (6.8) 

\f\=m 

6.2 Interior regularity of solutions 

Lemma 6.1 Let £1 C W 1 be a bounded Lipschitz domain. Pick two functions H,Z G Co°(f2) such 
that HZ = H, and assume that 

Y [A^BMOCfi) - 6 ( 6 - 9 ) 

\a\ = \P\=m 

where 

S < (6.10) 
pp' 

with a sufficiently small constant c(m, n, k) > 0. 

IfU e W™(U, loc) for a certain q < p and AU € Wp m (n, loc), then U G W™(Q, loc) and 

V~LU\\w™{n) < C (\\HA(-,D)U\\ w - m ^ + \\ZU\\ W rn^). (6.11) 

Proof. We start with a trick applied in [10J under slightly different circumstances. We shall use 
the notation Ay for the operator A(Y, Dx), where Y G and the notation $y for a fundamental 
solution of Ay in M. n . Then, with star denoting the convolution product, 

HU + 3>y * (A - A Y )(HU) = <$>y*{HAU) + <5>y*{[A } H]{ZU)) (6.12) 
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and, consequently, for each multi-index 7, [ — ^ | = m, 



Dt(HU) + Yl D a+ ^r * - A a p{Y))D^HU)) 

\a\ = \f3\=m 

= D^$y * {HAU) + L> 7 $y * {[A,H]{ZU)). 
Writing this equation at the point Y and using (|3,32|) . we obtain 



(6.13) 



(l-Cpp's) \\D^nu)\\ Lp{Q) 

|7|=m 

< C(p,K)(\\HAU\\ Wp - m{n) + \\[A,n)(ZU)\\ w - m{n) ). 
Let p' < n. We have for every V € W™{VL) 

[ ([A,H](ZU),V)dX = \A(HZU,V) - A(ZU,HV)\ 
Jn 



By Sobolev's theorem 



and 



< cCll^ll^m-l^JVll^m^) + \\ZU\\ W rn n (^HVH^m-l ^ ) 

1 n+p P n 



I^WH^m-l^ < c\\ZU\\ W m (fi) 

n-\-p 



I VI 



vi/^r 1 («) - c iiHiw™(n)- 

p'n v J P 



Therefore, 



f ([A,H](ZU),V)dX < C \\ZU\\ W m pn m \\V\\ W rn m 
J Q n+p 

which is equivalent to the inequality 



\\[A,H](ZU)\\ w - m{n) <c\\ZU\\ W rn (n) . 



(6.14) 



(6.15) 



(6.16) 



(6.17) 



(6.18) 



(6.19) 



In the case p' > n, the same argument leads to a similar inequality, where pn/(n + p) is replaced 
by 1 + e with an arbitrary e > for p' > n and e = for p' = n. 

Now, (j6,ll|) follows from (|6.14|) if p' > n and p' < n, q > pn/(n + p). In the remaining case the 
goal is achieved by iterating this argument finitely many times. □ 

Corollary 6.2 Let p > 2 and suppose that (fO|) and (KHUj) hold. If U € W^{n,loc) and AU G 
W~ m (n, loc), then U € Wp m (ft, loc) and 



\KU\\w™(p.) < c (\\Z A(-, D)U\\ w - m ^ + \\ZU\\ w rn-i^). 



(6.20) 
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Proof. Let Zq denote a real- valued function in Cq°(0) such that HZq = ri and ZqZ = Zq. By 

firm 

WKUWwrn^ <C(\\HA(-,D)U\\ w -rn^ + \\Z U\\ W rn(ty) (6.21) 

and it follows from (|f).8|) that 

\\Z U\\ 2 wr(n) <ck- x $IA{ZqU,Z q U). (6.22) 

Furthermore, 

\A(Z U,Z U) - A(U,Z$U)\ < cK- l \\ZU\\ wr i {n) WZdAWw^ny (6.23) 

Hence 

ll-^oWH^-m^ < ck 1 (\\Z AU\\ w - m ^ \\Zq W||vi/ 2 m (Q) + K 1 \\2M\\w£ l ~ 1 (p) \\ 2 o^\\w^(n)) (6.24) 
and, therefore, 

\\Z U\\ W rn^ < cn^dlZ AU\\ w -m^ + ti^WZUW^-i^). (6.25) 
Combining this inequality with ()6.21j) we arrive at (j6.20f) . □ 

6.3 Invertibility of A : V™' a (Q) — ► V~ m ' a (Q) 
Recall the infinitesimal mean oscillations as defined in (|1 . 10|> . 

Theorem 6.3 Let 1 < p < oo, < s < 1, and a = 1 — s — 1/p. Furthermore, let f2 be a bounded 
Lipschitz domain in W 2 . Suppose that the differential operator A is as in §6.1 and that, in addition, 

^ i^ap}osc(n) + {>}osc(cK7) < 8, (6.26) 
\a\=\/3\=m 

where 

PP+Jij^)^^ («■») 

for a sufficiently small constant c > independent of p and s. In the case m = 1 the factor 
8/s(l — s) in (|6.27|) can be replaced by 5. 
Then the operator 

A : V™' a (ft) — ► V p - m ' a (n) (6.28) 

is an isomorphism. 

Proof. We shall proceed in a series a steps starting with 

(i) The construction of the auxiliary domain G and operator C 
Let e be small enough so that 
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\a\ = \(3\=m 

for all balls in {B r }a with radii r < e and 



V / / \A a p(X)-A Q p(Y)\dXdY<25 (6.29) 
JB r nnJB r nn 



^Br 



v(X) - u(Y) daxdoY < 25 (6.30) 



/B r ndnJB r ndtt 

for all balls in {-B r }an with radii r < e. 

We fix a ball £? £ in {-B £ }on and assume without loss of generality that, in a suitable system of 
Cartesian coordinates, 

nnB £ = {X = {X', X n ) G B e : X n > (p(X')} (6.31) 

for some Lipschitz function tp : R n_1 — > R. Consider now the unique cube Q(s) (relative to this 
system of coordinates) which is inscribed in B £ and denote its projection onto R" --1 by Q'(s). Since 
Vty? = —v'/v n , it follows from (|6.3()j) that 



JBL JB 



Vip(X') - Vtp(Y') dX'dY 1 < c(n) 5, (6.32) 



r r 



where B' r = B r nK n_1 , r < e. Let us retain the notation p for the mirror extension of the function 
ip from Q'(e) onto R n_1 . 

We extend A a p from Q(e) n Q onto Q(e)\Q by setting 

A a p{X) := A af3 (X', -X n + 2<p{X')), X e Q(e)\n, (6.33) 

and we shall use the notation 2t a « for the periodic extension of A a n from Q{e) onto R™. 

Consistent with the earlier discussion in §5, we shall denote the special Lipschitz domain {X = 
(X',X n ) : X' £ W 1 ^ 1 , X n > (f(X')} by G. One can easily see that, owing to 2en _1//2 -periodicity 
of (p and A a p, 

^2 [^Ibmo(G) + [V^BMOOR™- 1 ) < c{n)5. (6.34) 

\a\ = \(3\=m 

Now, with the operator A(X,Dx) in $7, we associate an auxiliary operator C(X,Dx) in G given 
by (ET231) . 

(ii) Uniqueness. 

Assuming that IA £ lp m ' a (f2) satisfies CIA = in 17, we shall show that U € V^ m '°(0). This will 
imply that IA = which proves the injectivity of the operator (|6.28|) . 

To this end, pick a function Ti € Cq > (Q(e)) and write CiTLlA) = [£, Also, fix a small 

> and select a smooth function A on Ri_, which is identically 1 on [0, 1] and which vanishes 
identically on (2,oo). Then by (ii) in Lemma 15. 11 

C{HU) - [C, H] {A( Plcg /6)lA) € V 2 - m '°(G) n V~ m ' a (G). (6.35) 
Note that the operator 

[C, H}p-l : V p m > a (G) — > y p —> a (G) (6.36) 
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is bounded and that the norm of the multiplier p Teg A(p Teg /6) in V™' a (G) is 0(9). Moreover, the 
same is true for p = 2 and a = 0. 

The inclusion (j6.35j) can be written in the form 

C(HU) + M{ZU) G V- m ' a {G) n V~ m ' (G), (6.37) 
where Z G Cg°(M n ), ZTL = Ti and .M is a linear operator mapping 

V™' a (G) y- m ' a (G) and F 2 m '°(G) y 2 ~ m '°(G) (6.38) 

with both norms of order 0(6). 

Select a finite covering of f2 by cubes Qj(e) and let {TCj} be a smooth partition of unity 
subordinate to {Qj(e)}. Also, let Zj G C^°(Qj(e)) be such that Wj-Hj = TCj. By Gj we denote 
the special Lipschitz domain generated by the cube Qj(s) as in part (i) of the present proof. The 
corresponding operators C and Ai will be denoted by Cj and Aij, respectively. It follows from 
(IOT1) that 

HjU + XX -1 m j z i z k)(n k u) g y p m ' a (0) n v™ fi (n). (6.39) 

Taking into account that the norms of the matrix operator Cj Mj Zj Z k in the spaces V™' a (£l) and 
KT' U (Q) are 0(6), we may take 6 > small enoug h and obtain HjU G F 2 m,u (ft), i.e. U G F 2 m ' u (ft). 
Therefore, C : V^ m,a (0) — > Vp~ m ' a (Q) is injective. 

(iii) ^4 priori estimate. 

Let p > 2 and assume that W G V^ m,a (f2). Referring to Corollary 16.21 and arguing as in part (ii) of 
the present proof, we arrive at the equation 

Hj U + Yl^f M i Z i Zk ^ Hk U ) = ( 6 - 4 °) 

k 

whose right-hand side satisfies 

ll^llv?*-(n) < c(\\AU\\ vfm .a {n) + HWH^-i^j), (6.41) 

for some domain to with ZJ C £1. Since the V^' a (S7)-norm of the sum in ()6.4U|) does not exceed 
C0||W||y™."(Q), we obtain the estimate 

N| V ".»(n) <c(\\AU\\ v - m , a{n] + WUW^-x^). (6.42) 

(iv) End of proof . 

Let p > 2. The range of the operator A : V™' a (Q) — > V^~ m ' a (0) is closed by (j6,4()j) and the 
compactness of the restriction operator: V^ m ' a (f2) — > 14 7 2 m_1 (a;). Since the coefficients of the adjoint 
operator £* satisfy the same conditions as those of C, the operator C* : V^ l ' a (Q) — > V^ m ~ a (SX) is 
injective. Therefore, we conclude that C : 14, m ' a (f2) — » 14r m ' _a (f2) is surjective. Being also injective, 
£ is isomorphic if p > 2. Hence C* is isomorphic for p' < 2. This means that £ is isomorphic for 
p < 2. The result follows. □ 
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6.4 Traces and extensions 



Let C R n be a bounded Lipschitz domain and, for m 6 N as well as 1 < p < oo and — 1/p < a < 
1 — consider a new space, V7™' a (0), consisting of functions U £ L p (Q,loc) with the property 
that p a D a U £ Lp(£l) for all multi-indices a with |a| = m. We equip Wp ' a (0) with the norm 

\\U\\w™- a (n) '■= E \\ Da U\\L P (n, P (X)<*pdX) + \\U\\l p (w), (6-43) 

|a|=m 

where oj is an open non-empty domain, U C f2. An equivalent norm is given by the expression in 
()1.6|) . We omit the standard proof of the fact that 

C°°(Tt) ^ W™' a (ft) densely. (6.44) 

Recall that for p £ (1, oo) and s £ (0,1) the Besov space Bp(dQ) is then defined via the 
requirement l|1.7)) . If we introduce the L p -modulus of continuity 

w P (/,t):=(y / |/(X)-/(r)|fd ( 7 X day) 1/P , (6.45) 



|X-V|<t 
X,Y€BCl 



then 



B-m - 11/11 w + ( L ^S dt ) ^ > ( 6 - 46 ) 





uniformly for / € B p (d$l). 

The nature of our problem requires that we work with Besov spaces (defined on Lipschitz 
boundaries) which exhibit a higher order of smoothness. In accordance with [25], we now make 
the following definition. 

Definition 6.1 For p £ (l,oo), m £ N and s £ (0,1), define the (higher order) Besov space 
B p n ~ 1+S (317) as the collection of all finite families f = {/a}| a |<m-i of functions defined on dQ with 
the following property. For each multi-index a of length < m — 1 let 

R a (X,Y):=f a (X)- Yl J l fa+p{Y)(X-Yf, X,YedQ, (6.47) 

\/3\<m-l-\a\ 

and consider the L p -modulus of continuity 



r a (t):=lll \R a (X,Y)\Pda x da Y ) 1/P . (6.48) 



\X-Y\<t 
X,YedCt 



Then 



I roo r /£\p \l/p 

{ an)-= E ll/alli»(an) + E (/ ^ m -i+s-\ a \)+n-i dt ) < °°- ( 6 - 49 ) 



\a\<m—l \ot\<m— 1 
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For further reference we note here that for each fixed k > 0, an equivalent norm is obtained 
by replacing r a (t) by r a (Kt) in Q6.49[l . Also, when m = 1, the above definition agrees with (jl.7|) . 
thanks to fHH). 

A few notational conventions which make the exposition more transparent are as follows. Given 
a family of functions {f a }\a\<m-i on 9fl and X £ $7, Y, Z € <9S7, set 

P Q (X,Y):= ljfa+p(Y)(X-Yy, Va:[a|<m-1, 

|/?|<m-l-|«| ^ ( 6 _ 50 ) 

P(X,Y) :=P (CW)) (A,Y), 

so that 

R a (Y,Z) = f a {Y)-P a {Y,Z), Va: \a\<m-l, (6.51) 
and the following elementary identities hold for each multi- index a of length < m — 1: 



D? x P a (X,Y) = P a+f3 (X,Y), \p\<m-l 



a 



P a (X,Y)-P a (X,Z) = Yl R a+ p(Y,Z) {X p\ 

\P\<m-l-\a\ 



(6.52) 



See, e.g., p. 177 in [47J for the last formula. 

Lemma 6.4 For each 1 < p < oo, —l/p<a< 1 — 1/ 'p and s = 1 — a — 1/p, the trace operator 

Tr : W£> a (Q,) — ► B s p (dn) (6.53) 

is well-defined, linear, bounded, onto and has Vp' a (Q) as its null-space. Furthermore, there exists 
a linear, continuous mapping 

£ : B s p {8tt) — > W^ a (n), (6.54) 

called extension operator, such that Tr o £ = I (i.e., the operator \6.53\) has a bounded, linear 
right-inverse). 

Proof. By a standard argument involving a smooth partition of unity it suffices to deal with the 
case when Q is the domain lying above the graph of a Lipschitz function tp : R n_1 — ► R. Composing 
with the bi-Lipschitz homeomorphism R™ 3 (X',X n ) \— > (X' ,(p(X') + X n ) € further reduces 
matters to the case when $7 = R™ , in which situation the claims in the lemma have been proved in 

pr*T| □ 

We need to establish an analogue of Lemma l6.4l for higher smoothness spaces. While for = R™ 
this has been done by Uspenskii in p9| . the flattening argument used in Lemma 16.41 is no longer 
effective in this context. Let us also mention here that a result similar in spirit, valid for any 
Lipschitz domain but with B m - 1+s+1 / p (n) in place of W™' a (n) (cf. (|TT7t|1 for the relationship 
between these spaces) has been proved by A. Jonsson and H. Wallin in |26j (in fact, in this latter 
context, these authors have dealt with much more general sets than Lipschitz domains). The result 
which serves our purposes is as follows. 
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Proposition 6.5 Let 1 < p < oo, — 1/p < a < 1 — 1/p, s = 1 — a — 1/p G (0, 1) and m G N. Define 
the higher order trace operator 

tr m _! : W p m - a (0) — > (6.55) 

by setting 

tr m _i W := (fH Tr [D a ) , (6.56) 

L J a|<m— 1 

where the traces in the right-hand side are taken in the sense of Lemma \b\J\ Then \n\55\ )- Ji(i.5n]) is a 
a well-defined, linear, bounded operator, which is onto and has V™' a (£i) as its null-space. Moreover, 
it has a bounded, linear right-inverse, i.e. there exists a linear, continuous operator 

£ : B™~ 1+S (dn) — ► W™> a (n) (6.57) 

such that 



f = {fa} H <m-i^B m - 1+s (dn)^i\ a \Tr[D a (£f)} = f a , Va:|a|<m-1. (6.58) 

In order to facilitate the exposition, we isolate a couple of preliminary results prior to the proof 
of Proposition 16.51 

Lemma 6.6 Assume that ip : W 1 - 1 is a Lipschitz function and define $ : R n 1 — > d£l <^-> R n 

by setting <&(X') := (X' ,(p(X')) at each X' G M n_1 . Define the Lipschitz domain Cl as {X = 
(X',X n ) G R n : X n > if(X')} and, for some fixed m G N ; p G (l,oo) and s G (0,1) consider a 
system of functions f a G Bp(d£l), a G Ng, \a\ < rn — 1, with the property that 

F) n 

— [/ a (*(X'))] = ^/ Q+ej ($(X'))a fe <I>,(X'), 1 < fc < n- 1, (6.59) 

fc 3=1 

for each multi-index a of length < m — 2, where {ej}j is the canonical orthonormal basis in R n . 
Finally, for each I G {1, ...,m — 1} introduce Ai : = {(ii, —,ti) : < t\ < ■ ■ ■ < t\ < 1} ; and define 
R a (X, Y) as in j6'.^7| j. Then if a is an arbitrary multi-index of length < m — 2 and r := m — 1 — \ot\, 
the following identity holds: 

R a (<S>(X'), $(¥')) 

E {/ [/a+e, 1+ ... +e3r ($(F , +tr(^ , -^)))-/a +e , 1+ ... + e, r ($(^)) 

(n,...,j r )e{i,...,n} r jA r 

r 

x JJ V^ jk (Y' + t k (X'-Y'))-(X'-Y')dt r ---dt 1 j, X' , Y' G R n ~ x . (6.60) 

k=l 

Proof. We shall show that for any system of functions {f a }\ a \<m-i which satisfies (|6.59[1 . any 
multi-index a G Ng with \a\ <m — 2 and any I G N with I < r := m — 1 — \a\, there holds 
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\p\<i 



{ji,...Jl)e{l,...,n}' 



fa + e J1+ -+e n + t t (X' - Y'))) - fa+e^-.+e^ (*(^)) 



J] V$ jk (Y' + t k (X' - Y')) ■ (X' - Y') dt r -- dtxj. 



(6.61) 



k=l 



Clearly, ()6.6Uj) follows from H6.47() and (|6,61|) by taking I := r. 

In order to justify (|6.61|) we proceed by induction on I. Concretely, when I = 1 we may write, 
based on (|6.59|) and the Fundamental Theorem of Calculus, 



UHx')) - uhy')) - ^/ a+e .($(y0)(^(^0 - <^CO) 

3=1 



1 d_ r 
dt 



U*(Y' + t(X'-Y'))) 



dt 



3=1 



d r 



n , /-i r 

/ [f a+ej (Y' + t(X> - Y')) - f a+e M Y ')) 



3=1 



X 



(6.62) 



as wanted. To prove the version of Q6.61JI when / is replaced by I + 1 we split the sum in the 
left-hand side of (|6.61|) . written for I + 1 in place of I, according to whether |/3|<Zor|/3|=Z + l 
and denote the expressions created in this fashion by S\ and S2, respectively. Next, based on (|6.59|) 
and the Fundamental Theorem of Calculus, we write 



fa+e jl+ ...+e jl + t t (X' - Y'))) - f a+ e n+ -+e n (*(Y')) (6.63) 

= E / U + e jl+ - + e jl+ e i mY' + U +l (X' - Y')))V^(Y' + t l+1 (X' - Y 1 )) ■ (X' - Y') dt l+1 
i=l Jo 

and use the induction hypothesis to conclude that 



Si 



E \L f a+ e n+ - + e n+i mY' + t l+1 (X'-Y'))) 

H V$ jk (Y' + t k (X' - Y')) ■ (X' - Y') dt l+1 ■■■dt i y 



Ui,-,3l+l)e{l,...,n} l+1 1+1 



(6.64) 



fe=i 
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Thus, if 

Fj(t) :=<5> 3 {Y' + t{X' -¥'))-$> 3 {Y>), 1 < j < n, (6.65) 
we may express £1 in the form 



Sl= £ { / [/ Q+e , 1+ ... + e j!+1 (*(F / + t l+1 (X' -y'))) - /«+e il+ ...+e il+1 (*(y'))' 

(ji,...,j<+i)e{i,...,n} i+1 ■ 7Ai+1 
i+1 



X II v^ fc (y' + t fc (x , -y / ))-(x / -y , )^+i---dti} (6.66) 

fc=l 

+ Yl / a+e , 1+ ... +e , !+1 ($(y'))/ n^(t*)Ai+i-cfti. 

(ii,-,i!+i)e{i > ...,n}i+i ,/A *+ 1 fc=i 

Note that the first double sum above corresponds precisely to the expression in the right-hand 
side of (|6.61j) written with I replaced by I + 1. Our proof of (|6,61|) by induction is therefore complete 
as soon as we show that for each multi-index (3 of length I + 1, 

E / f[F; k (t k )dt l+1 ---dt 1 = ±(*(X')-$(Y')f. (6.67) 

1 )e{i,...,n}i+i JAl + 1 k=l P ' 



In turn, this is going to be a consequence of a general identity, to the effect that 



to fti J__ [ 



E L L L I[Fj k (tk)dtr--dt 1 = -F(to) (3 , (6.68) 



(i 1 ,... Ji )S{l,...,n}' " - 70 J0 fc=l 



for any Lipschitz function i* 1 = (Fx, F n ) : [0, 1] — > C n with ^(0) = 0, any point to £ [0, 1] any 
/ € N and any multi-index f3 of length I. Of course, the case most relevant for our purposes is 
when the -Fj's are as in ()6.65() . to = 1 an d when I is replaced by I + 1, but the above formulation 
is best suited for proving (|6,68|) via induction on /. Indeed, the case I = 1 is immediate from the 
Fundamental Theorem of Calculus and to pass from / to I + 1 it suffices to show that the two sides 
of (|6.68|) have the same derivative with respect to to- The important observation in carrying out 
the latter step is that the derivative of the left-hand side of (|6.68f) with respect to to is an expression 
to which the current induction hypothesis is readily applicable. This justifies (|6.68|) and completes 
the proof of (HOIUl) . □ 

Corollary 6.7 Under the assumptions of Lemma \6. 61 for each multi-index of length < m — 2 the 
following estimate holds 

/ , K1 ;-H)^i dt) ^ C E H/tIIw (6-69) 

7|=m— 1 

where the constant C depends only onn, p, s and \\^(p\\L ao (R n - 1 )- 
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Proof. The identity (|6.6Uj) gives 



\R a ($(X'),$(Y'))\ (6.70) 

< c\x' - yp-H«l Y I \M$( Y ' + T ( x ' ~ - MHY'))\dr 

, " Jo 



|7|=m— 1 



for each X',Y' E W 1 , where the constant C depends only on n and HV^Hl^ which, in turn, is 
controlled in terms of HV^H^. Given an arbitrary t > we now integrate the p-th power of both 
sides in (l6~m for X',Y' G dW 1 ' 1 subject to \$(X') - <$>{Y')\ < t. Using Fubini's Theorem and 
making the change of variables Z' := Y 1 + t{X' — Y') we obtain, after noticing that \Z' — Y'\ < rt, 

r a {tf < CtP^- 1 -^ Y J J J 1 ^{^{Y' + t(X' -Y')))- f^iY'WdrdX'dY' 

| 7 |=m-l x , t y, eRn -l ° 
\X'-Y'\<ct 

< ctP^- 1 -^ Y f o J J \fi($( z '))- MH Y 'W dZ'dY'dr 



7 =m— 1 



\Z' -Y'\<crt 



< Ct P(m-l +S -\a\)+n-l y f a; p (/ 7 , crtf 
Consequently, 



7|=m— 1 ' 



f r°W dt < C V !°° f 1 ^ CTt)P drdt 

J Q tp(m-l+s-\a\)+n-l ~ J Q J Q T n-ltf>r < 



'S+n—1 



7|=m— 1 

— ^ ^ > ( j n r ps+n-l ( j n T l-sp ^ T 



|7|=m— 1 







[7[=m— 1 ' 







after making the change of variables r := art in the second step. With this in hand, the estimate 
(KTM1) follows by virtue of (IPfil) . □ 

After this preamble, we are in a position to present the 

Proof of Proposition 16.51 We divide the proof into a series of steps, starting with 

Step I: The well-definiteness of trace. Let U be an arbitrary function in W p m ' a {Q) and set 

f a := TV [D a U], M a : \a\ <m - I. (6.73) 
It follows from Lemma 16.41 that these trace functions are well-defined and, in fact, 
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^ \\fa\\B°(dn) < C\\U\\ w ™,a( n y (6.74) 

\a\<m—l 

In order to prove that / := {fa}\ a \<m-l belongs to £>™ _1+s (9ft), let R a {X,Y) and r a (t) be as 
in l|6.47)) - (|6.48j) . Our goal is to show that for every multi- index a with \a\ < m — 1, 

' Q y(m -i + ° -|a| )+ n-i dt ) < C WU\\wr«(ny (6-75) 
To this end, we first observe that if \a\ = m — 1 then the expression in the left-hand side of 
(lfi~75|) is majorized by C(7 °° u p (f a , tf /t ps+n - 1 dt) P which, by (IQHj) and (jEHJ), is indeed < 
C 1 1 ^Hiy™' a (r2)- To treat the case when \a\ < m — 1 we assume that is locally represented as {X : 
X n > f(X')} for some Lipschitz function cp : M n_1 — > M and, as before, set := (X', (p(X')), 

X' G Then ()6.59|) holds, thanks to (|6.73j) . for every multi-index a of length < m — 2. 

Consequently, Corollary 1(1 71 applies and, in concert with (|fi.74j) . yields (j6.75j) . This proves that the 
operator (|6.55|) - (|6,56|) is well-defined and bounded. 

Step II: The extension operator. We introduce a co-boundary operator £ which acts on / = 
{fa}\ a \<m-i G B™- l+s (dQ) according to 

(£f)(X)= f fC(X,Y)P(X,Y)day, Ie!l, (6.76) 
Jan 

where P(X,Y) is the polynomial associated with / as in (|6.50[) . The integral kernel /C is assumed 
to satisfy 

f K(X, Y) day = 1 for all X G ft, (6.77) 
Jan 



'an 

\l—n—\a 



\D%lC(X,Y)\<c aP {Xy- n -W, VlGfi, VYeSfi, (6.78) 
where a is an arbitrary multi-index, and 

£pT,Y) = if \X-Y\>2p{X). (6.79) 
One can take, for instance, the kernel 

JC(X, Y) := r, ( X ~^ ) ( [ V ( X ~ ( l, ) daz) ' , (6-80) 
\xp Teg {X)J \J m \xp ies {X)J J 

where rj G C^°(i?2), ^ = 1 on B±, n > and x is a positive constant depending on the Lipschitz 
constant of 9ft. Here, as before, p reg (X) stands for the regularized distance from X to 9ft. 
For each X G ft and Z G 9ft and for every multi- index 7 with (7! = m we then obtain 

D^£f(X) = 7LI D a x K(X,Y){P p (X,Y)-Pp(X,Z))da Y . (6.81) 

H>i 

If for a fixed /i > 1 and for each X G ft and t > we set 

T t := {Y G 9ft : \X - Y\ < fit} (6.82) 
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we may then estimate 



\W£f(X)\* < C P(Xr pH -f \Pp(X,Y)-P p (X,Z)\ p da Y 

H>1 

< C E P{X)- pH j \R s+ p(Y,Z)\ P \X-Y\ pl5l da Y , 

H>1 ]yl ' '- 

< C p{X) pilThm) -f \R T (Y,Z)\ p da Y , (6.83) 

|r|<m-l r p( x > 

where we have used Holder's inequality and (jd52j) . Averaging the extreme terms in (|(183j) for Z in 
r p (v)> we arrive at 

\D~<£f{X)\ p < C p{XY^\' m ^ 2{ - n -^ j j \R T (Y,Z)\ p da Y da z . (6.84) 

|r|<m-l ^ r p(x)^r p(x) 

Consider now a Whitney decomposition of 0, into a family of dyadic cubes, {Qi}i<=i- In partic- 
ular, li : = diamQj ~ dist (Qi, d£l) uniformly for i £ Z. Thus, if X € Qi for some i € Ij := {i € X : 
h = 2~ j }, j £ Z, the estimate (|6~%4"|) yields 

|D 7 £/(X)|<C 2 -j(|T|-m)^ 2 2j(n-l)yy | J R r (y,Z)| p dayda Z ) 1/P (6.85) 

\Y-Z\<x 2~i 

for some x = x(d£l) > 1. In fact, by choosing the constant p in (|6.82() sufficiently close to 1, 
matters can be arranged so that the family {xQjjigx has finite overlap. Keeping this in mind and 
availing ourselves of the fact that p(X) ~ Zj uniformly for X € Qi, i € T, for each multi- index 7 of 
length m we may then estimate: 

\ D -ygf^X)\Pp(X) p ^- 1 dX 

1 

- C /ZY1 2~ M1 ~ s) ~ j [ \D~<£f(X) \ p dX 

- C /ZYj E 2^ m - 1+s -l r D + ^"- 1 ) j J \R T (Y,Z)\ p da Y da z 

jeZielj |r|<m-l Y.zednnxQi 

\Y~Z\<x 2-3 



- C Y1 Yj 2 jp(m - 1+s ~ lTl)+j{n ~ 1) \R T (Y,Z)\ p da Y da z 



< 



jeZ|i-|<m-l Y,zedn 

\Y-z\oc2-i 

c v r rAt)p dt 

|r|<m- 



/n ^p(m— 1+s— |r|)+n-l 
1 U 
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where in the last step we have used (|6.49j) . This proves that the operator (|6.54j) is well-defined and 
bounded. 

Step III: The right-invertibility property. We shall now show that the operator (|6,76|) is a right- 
inverse for the trace operator (j6.55|) . i.e., whenever / = {/ 7 }i 7 |< m _i G B™~ 1+s (dQ,), there holds 

/ 7 = »M Tr[D 7 £7] (6.87) 
for every multi-index 7 of length < m — 1. To this end, for I7I < m — 1 we write 

D-r£f(X)-Sj(X)= Y, I D%K{X,Y)(Pp{X,Y)-Pp{X,Z))do-Y, (6.88) 

a+/3 = 7 a -P- Jdil 

\a\>l 

where 

£ 1 f(X):=[ lC(X,Y)P 1 (X,Y)da Y , X G 0. (6.89) 
Jan 

Estimating the right-hand side in (|6.88|) in the same way as we did with the right-hand side of 
H6.81|) . we obtain 



/ \W£f{X)-& r f{X)\*p{X)-*'- 1 dX < C V / 

JdQ \r\<m-l J ° 



^ P dt 



£p(|7|+s-|i-|)+n-l 



< C\\f\\l rl+sm . (6.90) 



In a similar fashion, we check that 



/ \V(D^£f(X) - SjiXWpiXr^- 1 dX 
Jan 



1 an 

< 







|r|<m-l ' 

The two last inequalities imply D^£f — £f G V p 1,a (^l) and, therefore, 

Tr {D 1 £f - £ 7 f) = 0. (6.92) 

Going further, let us set 

Eg(X) := [ K{X,Y)g{Y)do Y , X G Q. (6.93) 
Jan 



A simpler version of the reasoning in Step II yields that E maps Bp(dQ) boundedly into Wp ,a (Q). 
Also, a standard argument based on the Poisson kernel-like behavior of IC(X, Y) shows that Tr Eg = 
g for each g G Bp(dft). 
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Based on the definition (|6.50|) and (|6.89j) we have 

\Syf(X) - Ef,{X)\P + p(Xy\V(£ 7 f(X) - Ef 7 (X))f 

<C Y, P( X ) m { \f 1+ p(Y)\ p da Y . (6.94) 

I/9|>1 

Consequently, for an arbitrary Whitney cube Qi we have 

/ \£j(X)-Ef 1 (XWp(X)-P s - l dX+ [ \V{£ 1 f{X)-Ef 1 {XWp{Xy-^- 1 dX 
JQi JBs 

<C Y, l i m ' S) [ \f 1+ p(Y)\ p da Y . (6.95) 

I/3|>1 

Summing over all Whitney cubes we find 

\\S-yf - Ef^\\ v i, a{n) < C ^ \\fa\\L p (dn) (6-96) 

|a|<m— 1 

which implies 

Tr(£j-Ef y ) = 0. (6.97) 
Finally, combining (|6.97|) . (|6.92|) . and TrEf^ = / 7 , we arrive at (|6.87|) . 

Step IV: The kernel of the trace. We now turn to the task of identifying the null-space of the trace 
operator (|6,55|) - H6.56|) . For each k € No we denote by Vk the collection of all vector- valued, complex 
coefficient polynomials of degree < k (and agree that Vk = whenever A; is a negative integer). 
The claim we make at this stage is that the null-space of the operator 

W™' a (ty B W ^ (ir [Dm]} € Bt(dn) (6.98) 

I J 7|=m— 1 

is given by 

V m - 2 + V p m ' a (tt). (6.99) 

The fact that the null-space of the trace operator (|6.55|) - ()6.56|) is V p m,a {Vt,) follows readily from this. 

That (|6.99f) is included in the null-space of the operator (|6.98|) is obvious. The opposite inclusion 
amounts to showing that if W G W™' a (£l) is such that Tr [Z? 7 W] = for all multi- indices 7 with 
I7I = m— 1, then there exists P m ~2 £ with the property that W — P m -2 £ V p m ' a (Q). To this 

end, we note that the case m = 1 is a consequence of (|5.19|) and consider next the case m = 2, i.e. 
when 

W G W^ a (n), Tr [VW] = on dO. (6.100) 

Assume that {VVj}j>i is a sequence of smooth in (even polynomial) vector-valued functions 
approximating W in Wp' a (Q). In particular, 
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Tr [VVWj] ^ in L p (dQ) as j -> oo. (6.101) 

If in a neighborhood of a point on d£l the domain is given by {X : X n > ip(X')} for some 
Lipschitz function </?, the following chain rule holds for the gradient of the function Wj : B' 9 I' h 
Wj(X',tp(X')), where B' is a (n — l)-dimensional ball: 



V Wj (X') = (v y ,W,(y>pn))| y/=x + (JLwj(X',Y n j)\ = V<p(X'). (6.102) 



.dY n 



Y n =ip(X>) 



Since the sequence {wj}j>i is bounded in L P (B') and Vu>j — > in L p (B'), it follows that there 
exists a subsequence such that Wj i — > const in L p (B') (see Theorem 1.1.12/2 in |31j). Hence, 
Tr W = Po = const on 50. In view of Tr [W — Po] = and Tr [VW] = 0, we may conclude that 
W - P G Vp 2 ' a (ft) by Hardy's inequality. 

The general case follows in an inductive fashion, by reasoning as before with D a W with \a\ = 
m — 2 in place of W. □ 

We now present a short proof of (jl.Hjl , based on Proposition 16.51 
Proposition 6.8 Assume that 1 < p < oo, s G (0, 1) and m G N. TTien 

\\f\\B™- 1+s (dn) ~ H/ttllB*(an) 5 (6.103) 

a|<m— 1 

uniformly for f = {f a }\ a 

<m— l G P^ 1 1 ^ s (9f2). As a consequence, holds. 

Proof. The left-pointing inequality in H6.103|) is implicit in (|6.69|) . As for the opposite one, let 
/ = {/«}| a |< m „ 1 G P™- 1+S (dft) and, with a := 1 - s - 1/p, consider U := £(f) G H/ p m ' a (0). Then 
Lemma 16.41 implies that, for each multi- index a of length < m — 1, the function f a = i ' a ' Tr [D a U] 
belongs to B p (dQ), plus a naturally accompanying norm estimate. This concludes the proof of 
(|6.103|) . Finally, the last claim in the proposition is a consequence of (|6,103j) . (|6.44j) and the fact 
that the operator l|6.55|) - l|6.56jl is onto. □ 

We include one more equivalent characterization of the space P™~ 1+S (<9f2), in the spirit of work 
in [lj, [42J, [50J. To state it, recall that {ej}j is the canonical orthonormal basis in W 1 . 

Proposition 6.9 Assume that 1 < p < oo, s G (0, 1) and m G N. Then 



f a G B°(dQ), Va : |a| < m- 1 
and 

(6.104) 

(Ujdf. — V k dj)f a = Vjf a+ek — fkfa+ej 

Va : \a\ < m — 2, V j, k G {1, ...,n}. 

Proof. The left-to-right implication is a consequence of (|6.103j) and of the fact that (|6.73|) holds 
for some U G W p ,a (fi) (cf. Proposition I6.5j) . As for the opposite implication, we proceed as in the 
proof of Proposition 16.51 and estimate (|6.48|) based on the identities Q6.59JI and knowledge that f a 
belongs to B p (dQ) for each a of length < m — 1. □ 



{/a}|a|<m-l G ^- 1+s (90) ^ < 
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We close this section with two remarks on the nature of the space B™ 1+s (<90). First, we claim 
that the assignment 

^ 1 L ' — » a! J 0<fc<m-l 

|cc|=fe 

is one-to-one. This is readily justified with the help of the identity 

^ = rH^|l+ £ Ep^mM-^-DP (6.106) 

|/3| = |a|-lj,fc=l ■ ?fc 

where d/drj^ := Vjd/dx}. — v^d/dxj and the p a ,p,j,kS are polynomial functions. Indeed, let 
/ 6 i?™~ 1+s ((?f2) be mapped to zero by the assignment 1)6.105(1 and consider U := £ (/) € W^"' a (fi). 
Then / a = i' Q ' Tr [D a Z^] on dQ. for each a with \a\ < m — 1 and, granted the current hypothe- 
ses, d k U/dv k = for k = 0,1,..., m — 1. Consequently, (|6.106|) and induction on |a| yield that 
Tr [D a U] = on d£l for each a with \a\ < m — 1. Thus, f a = for each a with |a| < m — 1, as 
desired. 

The elementary identity (|6.106j) can be proved by writing 



•|a| D a 



n 

3=1 



f— 

\dxi 



(6.107) 



n n 

n[s&(& 

3=1 k=l 

n OLj 



d 



8 



dxj dxk 



8 



k=l 




dx k 




II [E n( a /- /)! (X^fe, 
j=l «=o v 3 ' k=l 3 K 

n r 5 a J a \ / n / d 8 

n r? + E U r a 3 _ At (E ^ (^^: - ^ 

7=1 z=o v J ; fe=i 



5=^ 



and noticing that E[?=i v* J 8 a 3 / 8v a i = u a 8^ / 8v^ , whereas (^d/dxj — £jd/8xk)\ ^ =v = —8/drjk- 
Our second remark concerns the image of the mapping 1)6.105(1 in the case when d£l is sufficiently 
smooth. More precisely, assume that 80, G C m ~ 1,1 and, for < k < m — 1, the space B^ 1 ^ k+s (8U) 
is defined starting from £?^ _1— fe+s (R n_1 ) and then transporting this space to dfi via a smooth 
partition of unity argument and locally flattening the boundary (alternatively, B™~ 1 ~ k+S (80.) is 
the image of the trace operator acting from B™ 1 k+s+1 ^ p (W 1 )). We claim that 



dn G C™- 1 . 1 =>> the image of the mapping (KHUBl is ©J^B™- 1 "^'^). (6.108) 

Indeed, granted that dtt G C" 71 " 1 ' 1 , it follows from (l6~104l that f a G £™~ lHQ:|+s (dn) for each a 
with |a| < m — 1 and, hence, := i fc X)| a |=fc ST ^° e -B™ _1_fc+S (<9r2) for each fc G {0, m — 1}. 

Conversely, given a family {<?fc}o<fc<m-i G B™~ l ~ kJrS (8£l) , we claim that there exists 

/ = {f a }\ a \< m -i G B™- 1+S (dn) such that 9fc = i fc E| Q |=fc ST ^ /a for each fc G {0, ...,m - 1}. One 
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way to see this is to start with U G B™~ +s+ (f2) solution of A m U = in Q, d k U/dv k = g k G 
B™~ l ~ k+S (dfl) , < k < m — 1 (a system which satisfies the Shapiro-Lopatinskij condition) and 
then define the / Q 's as in l|6.73j) . 

6.5 Proof of the main result and further comments 

Theorem 11.11 is a particular case of the next theorem concerning the unique solvability of the 
Dirichlet problem in W™' a (£l). 

Theorem 6.10 Let all assumptions of Theorem \6.?\ be satisfied. Also let T G V p m ' a (ft). Then the 
Dirichlet problem 

' A(X,D X )U = T infl, 

- , = qu. on <9fL < k < m — 1, 

has a solution IA G W^"' a (0) if and only if U.y\) is satisfied. In this latter case, the solution is 
unique and satisfies 

\\U\\w™' a (n) < C ^2 WfaWs^dn) + C\\F\\ Vp -™^(n)- (6.110) 

|o|<m— 1 

Proof. It is clear from definitions that the operator 

A(X, D x ) : W; n ' a (Q) — » V~ m ' a (n) (6.111) 

is well-defined and bounded. Thus, granted that we seek solutions for (|6.109|) in the space W p ' a (Q), 
the membership of T to Vp~ rn ' a (Q), as well as the fact that the g^s satisfy p.9j) . are necessary 
conditions for the solvability of (|6.109|) . 

Conversely, let / = {f a }\a\ <m—i ^ -^p 1 ^~ >rS {d£l) be such that (|1.9|) holds and, with £ denoting 
the extension operator from Proposition 16. 5( seek a solution for (|6.109|) in the form IA = £(f) + W. 
where W G V™' a (Q) solves 

( A(X,D X )W = F-A(X,D x )(£(f)) in ft, 

{ (6-112) 
{ Tr [D~< W] = on dtt, V7 : |-y| < m - 1. 

Since the boundary conditions in (|6.112j) are automatically satisfied if W G V™' a {Q), the solvability 
of (|6.112|) is a direct consequence of Theorem 16.31 As for uniqueness, assume that IA G W™ ,a (Vt) 
solves (|6. 109(1 with T = and = 0, < k < m — 1. From the fact that (|6.105|) is one-to-one, we 
infer that Tr [D^IA] = on dtl for all 7 with | — y | < m — 1. Then, by Proposition 16. 5( IA G V™' a (£l) 
is a null-solution of A(X, D x ). In turn, Theorem 16.31 gives that U = 0, proving uniqueness for 
(|6.109|) . Finally, (|6.110|) is a consequence of the results in §6.4. □ 

We conclude this section with a couple of comments, the first of which regards the effect of the 
presence of lower order terms. More specifically, assume that 

A{X,D X )IA:= Yj D a {A a p(X)D l3 U), X £ Q, (6.113) 

0<\a\,\P\<m 
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where the top part of A(X, Dx) satisfies the hypotheses made in Theorem 1 1 . 1 1 and the lower order 
terms are bounded. Then the Dirichlet problem (|6.1U9j) is Fredholm solvable, of index zero, in the 
sense that a solution U G W™' a (Q) exists if and only if the data J 7 , {gk}o<k<m-i satisfy finitely 
many linear conditions, whose number matches the dimension of the space of null-solutions for 
(|6.1()9|) . Furthermore, the estimate 

\M\wr a m < c ( Yl WFWv-^M + WMhm) + \\uh P (n)) (6.114) 

\a\<m—l 

holds for any solution U G W™' a {Q) of ([B~TU9~j) . 

Indeed, the operator 

A : V™> a {Q) — » V p - m ' a {n) (6.115) 

is Fredholm with index zero, as can be seen by decomposing A = A + (A — A) where A := 
^2\ a \ = \p\ =m D a A a p D@ , and then invoking Theorem 16.31 Now, it can be shown that the problem 

()6.1(J9|) is solvable if and only if T— A(X, Dx)£f G Im„4, the image of the operator (|6.115|) . Thus, if 

T(F, {gk}o<k<m~i) ■= F-A(X,D x )£f, this membership entails (IF, {gk}o<k<m-i) G T^HmAj. 

Note that T maps the space of data onto Vp m,a {Q), hence the number of linearly independent 
compatibility conditions the data should satisfy is 

codimT -1 ^Im^ = codim (Im.4). (6.116) 

On the other hand, from by Proposition 16.51 and the fact that H6.105|) is one-to-one we infer that 
the space of null-solutions for (|6.1()9j) is precisely ker.4, the kernel of the operator (|6.115|) . Since, 
as already pointed out, this operator has index zero, it follows that the problem (|6. 109(1 has index 
zero. Finally, (|6.114|) follows from what we have proved so far via a standard reasoning as in [23J . 
Our last comment regards the statement of the Dirichlet problem (jl.2j) with data 

d k U/du k = g k G B™- l - k+s {dQ), fc = 0,l,...,m-l, (6.117) 

where £?™ _1_fc+s (<9f2) is defined here as the range of Tr acting from B™ 1 fc+s+1 / p (]j™), if [ s 
smooth (C ' will do) this problem is, certainly, well-posed. Let us illustrate some features of this 
particular formulation as the smoothness of d£l deteriorates. 

Suppose we are looking for the solution U G VF|(f2) of the Dirichlet problem for the biharmonic 
operator 

' A 2 U = in Q, 

< TrU = g 1 ondn, (6.118) 

(z/,Tr[W]) = g 2 on dtt. 

The simplest class of data (51,52) would be, of course, B^ /2 (dn) x Bl /2 (dn), where flj /2 (0n) and 
B^ 2 {dQ) are the spaces of traces on dfl for functions in W^^l) and W^i^l), respectively. However, 
this formulation has several serious drawbacks. 
The first one is that the mapping 

W 2 2 (0) B U -> (u, Tr [VW]) G -B 2 1/2 (5fi) (6.119) 



51 



is generally unbounded. In fact, by choosing U to be a linear function we see that the continuity 
of implies v E B^idVl) which is not necessarily the case for a Lipschitz domain, even for 

such a simple one as the square S = [0, l] 2 . 

The same problem fails to have a solution in the class in Wf(f2) when when (51,^2) is an 
arbitrary pair in B^ 2 (dQ) x b\ (cK2). Indeed, consider the problem (j(1118() for Q = S and the 
data g\ = and <?2 = 1. It is standard (see Theorem 7.2.4 in |29| and Sect. 7.1 in |3D|) that 
the main term of the asymptotics near the origin of any solution U in W^iS) is given in polar 
coordinates (r, u) by 

2r / , 7r \ 

ui — —) sin u> — 00 cos to I . (6.120) 



vr + 2 V 2 

Since this function does not belong to W^iS), there is no solution of problem ((6.118(1 in this space. 
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